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ABSTRACT 
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^ Aims. This paper describes a systematic search for high-order multiplicity among wide visual Pre-Main Sequence (PMS) binaries. 
■■^q- Methods. We conducted an Adaptive Optics survey of a sample of 58 PMS wide binaries from various star-forming regions, which include 52 
f"^ , T Tauri systems with mostly K- and M-type primaries, with the NIR instrument NACO at the VLT. 

Results. Of these 52 systems, 7 are found to be triple (2 new) and 7 quadruple (1 new). The new close companions are most likely physically 
bound based on their probability of chance projection and, for some of them, on their position on a color-color diagram. The corresponding 
degree of multiplicity among wide binaries (number of triples and quadruples divided by the number of systems) is 26.9 ±7.2% in the projected 
O ' separation range ~0'.'07-12", with the largest contribution from the Taurus-Auriga cloud. We also found that this degree of multiplicity is 
, twice in Taurus compared to Ophiuchus and Chamaeleon for which the same number of sources are present in our sample. Considering a 
(/} ■ restricted sample composed of systems at distance 140-190 pc, the degree of multiplicity is 26.8 ±8.1%, in the separation range 10/14 AU - 
' 1700/2300 AU (30 binaries, 5 triples, 6 quadruples). The observed frequency agrees with results from previous multiplicity surveys within 
the uncertainties, although a significant overabundance of quadruple systems compared to triple systems is apparent. Tentatively including 
*~J the spectroscopic pairs in our restricted sample and comparing the multiplicity fractions to those measured for solar-type main-sequence 
'Nj ' stars in the solar neighborhood leads to the conclusion that both the ratio of triples to binaries and the ratio of quadruples to triples seems 
^ to be in excess among young stars. Most of the current numerical simulations of multiple star formation, and especially smoothed particles 
' hydrodynamics simulations, over-predict the fraction of high-order multiplicity when compared to our results. The circumstellar properties 
around the individual components of our high-order multiple systems tend to favor mixed systems (i.e. systems including components of 
wTTS and cTTS type), which is in general agreement with previous studies of disks in binaries, with the exception of Taurus, where we find a 
preponderance of similar type of components among the multiples studied. 
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1. Introduction 

It has long been recognized that the formation of binary or mul- 
tiple stars is an efficient way to solve one of the most prominent 
problems in the theory of star formation, namely the "angu- 
lar momentum problem" (Mestel & Spitzer 1956] Larson [T972l 
Mouschovias l 19771 Simon et al l 19951 Bodenheimei l 19951 
Larson 2002 1. Yet, although we have a general understand- 
ing of the formation of single low-mass stars (e.g. Shu et 
al. fl987b . the theory of the origin of binary and multiple low- 
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mass stars is not quite settled (for a review see Tohline2002i. 
The favoured mechanism for the formation of binary and 
multiple stars is the collapse and fragmentation of molecular 
cloud cores, either fragmentation during dynamical collapse 
(e.g., Boss 1986 1 or from a quasistatic situation (i.e. disk frag- 
mentation e.g., Bonnell 1994 1 or filament fragmentation (e.g. 
Zinnecker Tl991l Bonnell & Bastien lT99"TV Prompt fragmenta- 
tion following the collision of two Jeans masses has also been 
suggested (Pringle 1991 1. In addition, further processes likely 
to occur during the formation of binary and multiple stars 
like circumbinary disk accretion, interaction with circumstellar 
disks and dynamical interactions are nowadays being investi- 
gated in detail through numerical simulations (Bate et al. 120031 
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Sterzik & Durisen l2UU51l2U031 Goodwin et al. EDCT . The scope 
of these studies is mainly to obtain statistical properties of the 
stellar systems resulting from these evolutionary calculations 
for direct comparison with observations; such as, relative fre- 
quency of single, binary and multiple stars, the properties of 
multiple stars and brown dwarfs, and the initial mass function. 

Several multiplicity surveys among T Tauri stars using 
high-angular resolution techniques have been carried out over 
the last decade, with the main outcome that binary stars are 
much more common in several of the young star-forming T- 
associations compared to the main-sequence stars in the solar- 
neighborhood (Leinert et al. lT993l Ghez et al. fT993l Ghez et 
al. U997al Duchene|l999i. In addition, it seems that a compa- 
rable or even higher frequency of binary and multiple stars 
is found at earlier evolutionary stages (e.g. Reipurth 2000, 
Duchene et al. 2006 and references therein). Paradoxically, al- 
though studies of high-order multiple star (N>2) systems might 
provide important constraints about the process of multiple star 
formation (Batten fT973l Fekel fT98l1 Sterzik & Tokovinin l2002l 
Tokovinin & Smekhov2002, see also Toko vinin 200 1 1, rela- 
tively little has been done observationally in a systematic way 
regarding higher-order multiplicity among Pre-Main Sequence 
(PMS) stars. An exception is Koresko (2002i who observed 
the close environment of 14 PMS binaries in southern star- 
forming regions with high-angular resolution techniques and 
found that half of them are in fact hierarchical triples or present 
resolved circumstellar structures consistent with marginally re- 
solved companions. In the last years, an increasing number of 
high-angular resolution observations reports the discovery of 
new companions in known multiple T Tauri systems. Examples 
are T Tau itself ( Koresk o lIUUOl Kohler et al. lSOOOat, VW Cha 
(Brandeker et al. |200"T1 . V773 Tau (Duchene et al.|2002), AS 
353 (Tokunaga et al.2004i. Some of these systems present 
close to edge-on circumstellar disks, more easily detectable in 
multiple systems and with current imaging techniques like AO. 
The only known cases of such disks are HK Tau B (Stapelfeldt 
et al. [T9^8l Kore sko[T998t . HV Tau C (Monin & Bouvie rEOOOl 
Stapelfeldt et al. l2003l, an d LkHa 263 C (Chauvin et al. l2002l 
Jayawardhana et al. 2002 1. Also, it has been noted that approx- 
imately half of the known T Tauri spectroscopic binaries are in 
fact hierarchical triples with an outlier companion (Melo 2003 
Sterzik et al. 2005 1. All the above suggests that with the current 
observational techniques we are able to reveal a large enough 
number of high-order multiple PMS stars to allow for mean- 
ingful statistical analysis of their occurrence and properties. 

We report here about a high-angular resolution survey 
which focuses on the N>2 multiplicity among T Tauri 
stars. We searched in a sample of 58 wide PMS bina- 
ries extracted from the list of Reipurth & Zinnecker ( 1993 
hereafter RZ93) to discover additional close companions 
using the Adaptive Optics (AO) near-infrared instrument 
NAOS/CONICA (NACO, Rousset et al. !2UU2l attached to the 
telescope UT4/Yepun of the Very Large Telescope (VLT). The 
rationale behind this effort was the following : if one consid- 
ers a wide 2 arcsec PMS binary (corresponding to ~ 300 AU 
at 150pc, the distance of the most prominent star forming re- 
gions), then a close (~ 0.1-0.2 arcsec) companion could still ex- 
ist without violating the hierarchical triple system stability cri- 



terion (Harrington[l975 Eggleton & Kiseleva n'9951 Mardling 
& Aarseth200TJl. 

This paper is organized as follows : in Sect.|2 we present 
the sample of wide PMS binaries and the AO near-infrared ob- 
servations obtained at VLT with NACO in Sect.[3] Candidate 
visual companions to these binaries detected in the images are 
presented in Sect.@]and physical companions are identified on 
both colors and statistical grounds. This allows us, in Sect. 14.71 
to derive the frequency of triples/quadruples systems among 
wide young binaries and to compare this quantity to previous 
surveys (Sect. l5.3l . to recent results of numerical simulations of 
multiple star formation (Sect. l5.4t . as well as to what is found 
among main-sequence multiples (Sect. l5.5l and 15. 6> . We also 
discuss the properties of these triples/quadruples systems in 
terms of dynamical stability (Sect. l5.8t and relative disk evo- 
lutions (Sect. l5.9l . Finally, we summarize our main results in 
Sect.© 



2. The sample 

We observed a sample composed of 58 wide PMS binaries ran- 
domly chosen from the 87 listed in RZ93 (see Table0. The 
separations range from 0f.'5 to 16'.'7, with a median of 2'.'4 and 
the bulk of the binaries having separations $ 6" (Fig^. It in- 
cludes sources from various nearby southern star-forming re- 
gion such as Chamaeleon, Gum Nebula, Lupus and Ophiuchus, 
while northern targets are mainly from Taurus-Auriga. Some 
objects from RZ93 and included in our sample are now known 
to be located closer and belong to the TW Hya and MBM12 
associations. Nearly all sources are T Tauri binary stars, hence 
composed of young, late-type (mainly K and M-type here) 
and low mass stars. Five sources (HD 76534, Herschel 4636, 
PHa30, Sz 120, Elias 2-49) are of earlier-type, intermediate 
mass Herbig Ae/Be binary stars. The sample we consider here 
is presumably strongly biased by selection effects which are 
difficult to quantify because they are the results of a com- 
plex combination of discovery biases. The sample is obviously 
not volume-limited since the original sample of RZ93 was 
mainly drawn from several Ha surveys which are magnitude- 
limited. This might have favoured the detection of multiple sys- 
tems which are brighter than single stars, or in our case, high- 
order multiples which are brighter than binaries (Opikl924i. 
However, we note that our sample is less biased towards bright 
systems in comparison to those of previous multiplicity surveys 
(see Sect. l5.3t . Our sample tends also to include mostly actively 
accreting objects, i.e. Classical T Tauri systems (CTTS). We 
should also point out that, since the survey of RZ93, sparse NIR 
high-angular resolution observations have either confirmed or 
discovered that ten wide binaries of the sample have additional 
physical companions. These systems are identified in the col- 
umn 'additional visual companions' in Tabled In addition to a 
likely triple system (Sz 30), two other systems were reported in 
RZ93 as possible triple systems (LkHff 336 and Sz 41). Despite 
all those known and/or suspected multiple systems, not all 
sources of our sample have been searched for additional com- 
panions (especially the faintest one) and our present AO survey 
uncovered additional new candidate companions. 
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Name 


HBC 


Other 


R.A. 


Decl. 


Cloud 


D 


V 


K 


RZ93 1 


add. 


SB 


SpTl 


SpT2 


Ref 






name 


[J2000.0] 




[pc] 






P 


PA 


comp." 








(SpT) 


Lktto 262/263 .. 


8+9 




02 56 08.4 +20 03 40 


MBM 12 


275 


14.6 


9.5 


4.1 


57 


X 




M0 


M0+M4 


1 


J 4872 






04 25 17.6 +26 17 51 


Tau-Aur 


142 


13.0 


8.6 


3.5 


231 


X 




K9 


Ml 


2 


FVTau 


386+387 




04 26 53.6 +26 06 55 


Tau-Aur 


142 


15.4 


7.4 


0.7 


289 


X 




K5+K6 


M3+M5 


1 


UX Tau 


43+42 




04 30 04.0 +18 13 49 


Tau-Aur 


142 


10.7 


7.6 


2.7 


270 


X 




K4 


M2+M3 


1 


DK Tail 


45 




04 30 44.3 +26 01 25 


Tau-Aur 


142 


12.6 


7.1 


2.4 


1 19 






K7/K9 


K7/M1 


3,2 


HK Tau 


48 




04 31 50.6 +24 24 18 


Tau-Aur 


142 


15.0 


8.6 


2.4 


1 72 






Ml 


M2 


4. 2 


LkHo- 266 


51+395 


V710Tau 


04 31 57.8 


+ 18 21 37 


Tau-Aur 


142 


14.6 


8.5 


3.2 


177 






M0.5 


M2 


3 


GG Tau 


54 




04 32 30.3 +17 31 41 


Tau-Aur 


142 


12.3 


7.4 


1.4 


185 


X 




K7+M0.5 


M5.5+M7.5 1 


UZ Tau 


52+53 




04 32 43.0 +25 52 32 


Tau-Aur 


142 


12.9 


7.4 


3.6 


275 


X 


1 


Ml 


IVIZ+IVIZ 


1 


HN Tau 


60+406 




04 33 39.3 +17 51 53 


Tau-Aur 


142 


13.7 


8.4 


3.2 


217 






K5 


M4 


j 


IT Tau 




Haro 6-26 


04 34 13.9 +26 11 42 


Tau-Aur 


142 


14.9 


7.9 


2.4 


224 






K3 


M4 




L1642-1 


413 


EWEri 


04 35 02.3 


-14 13 41 


LI 642 


140 


13.7 


7.7 


2.7 


349 






K7 




5 


RW Aur 


80+81 




05 07 49.6 +30 24 05 


Tau-Aur 


142 


10.3 


7.0 


1.4 


254 




2 


Kl: 


K5' 


4 1 


CO Ori 


84 




05 97 ^8 1 


+ 11 25 39 


Orion 


460 


10.6 


6.5 


2.0 


280 






F8:e 




4 


AROri 


155 




05 35 54.1 


-05 04 14 


Orion 


460 


13.9 


9.8 


2.0 


249 






K7/M0 




4 


LkHff 336 


190+516 




05 54 20.1 


+01 42 56 


LI 622 


460 


14.4 


9.2 


5.8 


95 






K7 


M0 5 




CGHff 5/6 






07 31 37.4 


-47 00 22 


Gum Neb. 


450 


14.2 


9.1 


11.1 


317 






K7 


K2/K5 


O 


)•[]« 14 


554 




08 08 33.8 


-36 08 10 


Gum Neb. 


450 


15.8 


10.3 


0.6 


164 






M2: 




4 7 


PHo- 30 




SPH42 


08 12 05.6 


-35 31 45 


Gum Neb. 


450 


15.1 


12.2 


0.7 


345 






B2e 




$ 


vBH 16 




vdBH 16 


08 27 39.0 


-5109 50 


Gum Neb. 


450 


15.6 


9.3 


4.5 


330 












PHiv 51 


561 




08 15 55.3 


-35 57 58 


Gum Neb. 


450 


15.9 


11.1 


0.6 


185 






K7/M0 




8 


HD 76534 




He 3-225 


08 55 08.7 


-43 28 00 


Vela R2 d 


830^ 


8.0 


7.8 


2.1 


303 






B2 


B^ 


14 


SX Cha 


564 


T3 


10 55 59.9 


-77 24 41 


Cha I 


160 


14.7 


8.7 


2.2 


288 






M0.5 




4 


Sz 15 




T19 


11 05 41.5 


-77 54 44 


Cha I 


160 


13.2 


10.6 


10.6 


28 






K5/K6 




7 


ESO Ha 281 






11 07 04.0 


-76 31 45 


Cha I 


160 


13.7° 


9.7 


1.7 


167 






M0.5 


M4.5 


10 


Sz 19 


245 


DI Cha 


1 1 07 20.7 


-77 38 07 


Cha I 


160 


10.9 


6.2 


4.6 


203 






G2 




4 


VV Cha 


573 


T27 


1 1 07 28.4 


-76 52 12 


Cha I 


160 


14.8 


9.5 


0.8 


14 






M1.5 


M3 


10 


VW Cha 


575 


T30/T31 


1108 01.8 


-77 42 29 


Cha I 


160 


12.6 


7.0 


16.7 


221 


X 


3 


K5/K7 


K7 


1 


Glass I 




Ced 1 1 1 IRS 4 


11 08 15.4 


-77 33 54 


Cha I 


160 


13.3 


6.9 


2.4 


105 






K3 


G5e 


11 


CoD -29°8887 .. 




TWA 2 


1 1 09 14.0 


-30 01 39 


TWHya 


50 


1 1.1 


6.7 


<0.8 


40 






M2e 


M2 


12 


Sz 30 




T39 


11 09 12.3 


-77 29 12 


Cha I 


160 


13.2 


9.0 


1.1 


22 






M0.5 


M2 


10 


Hen 3-600 




TWA 3 


11 10 28.9 


-37 32 05 


TWHya 


50 


12.1 


6.8 


1.4 


230 




4 


M3e 


M3.5 


12 


Sz41 


588 


T51 


11 12 24.5 


-76 37 06 


Cha I 


160 


11.6 


8.0 


1.9 


164 




5 


K3.5 




9 


CV Cha 


247+589 


T52 


11 12 27.8 


-76 44 22 


Cha I 


160 


1 1.0 


6.9 


11.4 


105 




6 


G8 


G: 


4 


Sz 48 






13 00 53.2 


-77 09 10 


Cha II 


178 


17.4 


9.5 


1.5 


232 






K7/M0 


M0 


10 


BK Cha 






13 07 09.3 


-77 30 24 


Cha II 


178 


15.2 


8.4 


0.8 


329 






K5/K7 


M0.5 


10 


Sz 60 






13 07 23.4 


-77 37 23 


Cha II 


178 


17.1 


9.5 


3.4 


280 






Ml 


M4 


13 


Sz 62 






13 09 50.7 


-77 57 24 


Cha II 


178 


15.6 


9.1 


1.1 


261 






M2 


M3.5 


10 


Herschel 4636 ... 




Hen 3-949 


13 57 44.1 


-39 58 45 


NGC 5367 


630 


9.7 


7.2 


3.7 


214 






B3 




14 


ESO Ho- 283 






15 00 29.6 


-63 09 46 


Circinus 


700 


15.5 


10.3 


2.1 


242 












Sz 65 


597 




15 39 27.7 


-34 46 17 


Lup I 


190 


12.7 


8.0 


6.4 


98 






MO 




15 


Sz 68 


248 


HT Lup 


15 45 12.9 


-34 17 31 


Lup I 


190 


10.4 


6.5 


2.6 


295 


X 




K2 


M6 


1 


HO Lup 


612 




16 07 00.6 


-39 02 19 


Lup III 


190 


13.0 


8.6 


1.5 


36 






K7/M0 


M2 


10 


Sz 101 






16 08 28.4 


-39 05 32 


Lup III 


190 


15.5 


9.4 


0.8 


305 






M2.5 


M3.5 


10 


Sz 108 


620 




16 08 42.7 


-39 06 18 


Lup III 


190 


13.1 


8.8 


4.2 


25 






Ml 




15 


Sz 120 






16 10 10.6 


-40 07 44 


Lup III 


1 90 


7. 1 


6.2 


2.6 


142 






B5 




15 


WSB 3 






16 18 49.5 


-26 32 53 


Oph 


160 


15.0 


9.3 


0.6 


162 






M3 




16 


WSB 11 






16 21 57.3 


-22 38 16 


Oph 


160 


18.5° 


10.1 


0.5 














WSB 20 


257 


Haro 1-4 


16 25 10.5 


-23 19 14 


Oph 


160 


13.4 


7.5 


<1.0 


23 






K6/K7 




4 


WSB 28 




ISO-Oph 27 


16 26 20.7 


-24 08 48 


Oph 


160 


14.7° 


9.5 


5.1 


358 






M3 


M7 


17 


SR 24 


262 


WSB 42 


16 26 58.8 


-24 45 37 


Oph 


160 


15.9 


7.1 


5.2 


348 


X 




K2: 


M0.5 


1 


Elias 2-30 




SR21 


16 27 10.2 


-24 19 16 


Oph 


160 


14.1 


6.7 


6.4 


175 






G2.5 


M4 


17 


WSB 46 


641+640 


ROXs 20 


1627 15.1 


-24 51 39 


Oph 


160 


14.9 


9.4 


10.3 


301 






M0/M3 


M1/M4 


18 


Haro l-14c 


644+267 




16 31 04.4 


-24 04 32 


Oph 


160 


12.3 


7.8 


12.9 


122 




7 


K3 


M0 


4 


ROXs 43 




NTTS 162819-2423 


16 31 20.1 


-24 30 05 


Oph 


160 


10.6 


6.7 


4.3 


13 


X 


8 


GO 


K3 


1 


Elias 2-49 




HD 150193 


16 40 17.9 


-23 53 45 


Oph 


160 


8.9 


5.5 


1.1 


227 






Al/A3Ve 




19 


HBC 652 


652 


V2507 Oph 


16 48 18.0 


-14 11 15 


L162 


160 


13.5 


7.5 


8.7 


322 






K4/K5 




4 


LkHff 346 


275 


THa 27-4 


17 11 03.9 


-27 22 57 


B59 


160 


14.3 


8.1 


5.0 


329 






K7 


M5 


1 



Note : HBC = catalog entry number in Herbig & Bell 1988 
a : R magnitude. 

b : Additional visual companion(s) discovered since RZ93 (see references in TableHlor in the text). 

c : Astrometr y of th e wide binary as quoted in RZ93 with p in arcsec and PA in degree (only one entry for the possible triples quoted in RZ93). 
d : see Herbsl 1975 and van den Ancker et al lT^SI 

Spectral type referen ces : (1) se e TableQQ (2) Duchen elT^I (3) Hartigan et al JT^I (4) Herbig & BelfTOfl (5) Sandel l et allp?S71 (6) Reipurth & Pe tterssortT^l (7) Kim et al|2005l 
(8) PetterssorHWp (9) Luhman 2004 (10) Brandner & Zinneckei 199^ (11) Feigelson & Kri s JOT)! (12) Torres et alJ2003l (13) Hughes & Hartiga ntT^l (14) Viera et al J2003l (15) 
Hughes et al ll9W (16) Meyer et al jl993l (17) Prato et al ]2TOl (18) Luhman & RiekeJTg?! (19) Gray & Corbally [T9981 

Distance references : Cham aeleon (1 : 160pc ± 15pc,II: 178 pc ± 18pc. Whittet et alffi?7l. Gum Nebula (450 pc, Kim et alJSTOH, Lupus (190 pc ± 27 pc, Wic hmann et al il998l . 
Ophiuchus (160pc, Whittet HTOI, Taur us-Auriga (142pc ± 14pc, Wichmann et al llM . TW Hya (50pc ± 12pc, Mamaiel TOBl . MBM 12 (275pc, Luhma nl2TjoTl . L1642, aka 
MBM20 (1 12-160 pc, Hearty et al l2TOI 140 pc adopted), NGC 5367, Cir cinus, L162 and B59 (values reported in RZ93 adopted). 

SB notes: (1) UZ Tau A : 19.1-days single- lined SB (SB1) (Mathi eu et a l J1996l . converted to double-lined SB (SB2) by Prato et al J2002l (2) RW Aur A: suspect ed SB 1 with P=2.77- 
days (Gahm et al lT939l Pet rov et al l^oTTTl see also Alencar et al l2TCl . (3) VW Cha A : possible SB2 (Melr CTBl . (4) TWA 3A : suspected SB2 (Torres et al lSpTBl . (5) Sz 41 A : 
probable SB (Reipurth et a U2002l with period of about 125 days (best candidate out of five radial-velocity variables). (6) CV Cha: possibly SB, variable radial-velocity (Reipurth et 
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Table 2. Observation Log. 



Name 


Obs. date 


Cam 


Int. time [s 


1 


WFS-Dic. 


seeing" 
[arcsec] 


FWHM 


mas] 


Airmass 


Strehl ratio [%] 


[Fell] 




Bry 


[Fell] 




Bry 




[Fell] H 2 


Bry 


LkHff 262/263 .. 


2002 Nov 14 


S27 


15x2.7 


15x1.5 


15x1.4 


VIS-VIS 


1.0 


129 


89 


88 


1.43 


4 


14 


12 


J 4872 


2002 Nov 13 


S13 


30x1.1 


30x1.6 


30x1.6 


VIS-VIS 


0.8 


72 


73 


72 


1.62 


18 


42 


40 


FVTau 6 


2002 Nov 13 


S27 


9x47 


9x19 


9x19 


IR-N20C80 


0.8 


121 


89 


81 


1 .59 


7 


19 


1 3 


UX Tau 


2002 Oct 22 


S13 


60x0 6 


60x0 6 

\.)\J AU.U 


60x0.6 


VIS-VIS 


1.1 


60 


68 


68 


1.40 


22 


38 


47 


DKTau 


2002 Nov 13 


S13 


30x1 .4 


30x 1 .0 


30x0.9 


IR-N20C80 


0.9 


64 


70 


68 


1 .59 


15 


29 


32 


HK Tau 


2003 Feb 17 


S13 


12x5.5 


12x2.9 


12x2.9 


IR-N20C80 


1.0 


SO 


73 


73 


1.75 


7 


21 


26 


LkHff 266 


2002 Oct 22 


S27 


105x0.3 


105x0.3 


105x0.3 


VIS-VIS 


1.7 


90 


82 


93 


1.44 


5 


12 


$ 


GG Tau 


2002 Oct 22 


S27 


60x 1 .0 


30x 1 .0 


30x1.0 


VIS-VIS 


1.7 


67 


70 


87 


1.49 


12 


31 


19 


UZ Tau 


2002 Nov 14 


S13 


45x0 6 


45x0 5 


45x0.5 


VIS-VIS 


1.0 


84 


nn 


77 


1.68 


8 


17 


27 


HN Tau 


2003 Feb 18 


S13 






15x1.5 


VIS-VIS 


0.6 






72 


1 .59 






29 


IT Tau 


2003 Feb 19 


S13 


15x3.8 


15x1.7 


15x1.7 


IR-N20C80 


1.1 


104 


71 


75 


1.84 


2 


17 


21 


L1642-1 


2003 Feb 20 


S13 


60x0.9 


60x0.6 


60x0.6 


VIS-VIS 


1.1 


78 


72 


77 


1 .50 


7 


23 


22 


RW Aur 


2002 Nov 17 


S13 


45x0.7 


45x0.5 


45x0.5 


IR-N20C80 


0.8 


66 


70 


70 


1.79 


15 


34 


36 


CO Ori 


2002 Nov 13 


S13 


30x2.0 


30x1.0 


30x2.0 


VIS-VIS 


1.0 


68 


71 


73 


1 .53 


12 


23 


28 


AR Ori 


2003 Feb 19 


S13 


3x20 


3x10 


3x10 


VIS-VIS 


1.1 


151 


122 


91 


1.21 


3 


8 


11 


LkHff 336 


2003 Feb 19 


S27 


6x10 


6x9.0 


6x8.5 


VIS-VIS 


1.4 


89 


87 


90 


1.20 


6 


12 


12 


CGHff 5/6 


2002 Dec 23 


S27 


36x1.3 


36x1.0 


36x1.0 


VIS-VIS 


0.7 


85 


84 


79 


1.08 


6 


13 


21 


Pita 14 


2002 Dec 23 


S27 


15x4.0 


15x2.0 


15x2.0 


VIS-VIS 


0.7 


86 


80 


78 


1.02 


8 


16 


24 


PHff 30 


2003 Jan 17 


S27 


3x28 


3x21 


3x20 


VIS-VIS 




104 


84 


83 


1.02 


9 


27 


32 


vBH 16 


2003 Jan 27 


S13 


15x5.5 


15x4.5 


15x4.5 


IR-N20C80 


0.9 


122 


112 


97 


1.15 


4 


10 


13 


PHa 51 


2003 Jan 23 


S13 


3x85 


3x31 


3x30 


VIS-VIS 


1.0 


135 


107 


88 


1.02 


3 


7 


11 


HD 76534 


2003 Jan 20 


S13 


60x0.6 


60x0.6 


60x0.6 


VIS-VIS 


1.6 


69 


80 


72 


1.18 


1 2 


1 8 


26 


SX Cha 


2003 Jan 22 


S13 


12x8.5 


12x2.5 


12x2.5 


IR-N20C80 


0.9 


93 


82 


76 


1 .66 


7 


J 


21 


Sz 15 


2003 Jan 22 


S13 


12x3.9 


12x4.0 


12x3.9 


VIS-VIS 


0.7 


85 


85 


83 


1.67 


7 


1 7 


16 


ESOHff281 


2003 Jan 22 


S13 


3x51 


3x34 


3x33 


VIS-VIS 


0.8 


99 


83 


87 


1.62 


5 


1 6 


16 


Sz 19* 


2003 Jan 20 


S13 


6x22 


6x16 


6x16 


VIS-VIS 


1.5 


69 


76 


76 


1 .66 


1 2 


22 


23 


VV Cha 


2003 Feb 19 


S13 


3x18 


3x8.0 


3x7.5 


IR-N20C80 


0.8 


173 


108 


100 


1 .66 


3 


9 


9 


VW Cha 


2003 Feb 20 


S27 


75x0.8 


75x0.3 


75x0.3 


IR-N20C80 


1.2 


77 


76 


72 


1.68 


1 1 


25 


33 


Glass I 


2003 Feb 20 


S13 


12x5.0 


12x3.5 


12x3.0 


IR-N20C80 


1.3 


88 


79 


82 


1.67 


5 


z 1 


24 


CoD -29°8887 .. 


2003 Jan 20 


S13 


90x0.3 


90x0.3 


90x0.4 


VIS-VIS 


1.1 


62 


69 


69 


1.01 


1 8 


a | 


34 


Sz 30' 


2003 Jan 20 


S13 






3x20 


IR-N20C80 


1.5 






184 


1.67 






7 


Hen 3-600 


2003 Feb 20 


S13 


60x0.4 


60x0.4 


60x0.4 


IR-N20C80 


1.5 


74 


91 


72 


1 .03 


10 


19 


26 


Sz 41 


2003 Jan 22 


S27 


75x0.8 


75x0.5 


75x0.5 


VIS-VIS 


0.9 


70 


74 


78 


1.64 


I ] 


27 


30 


CV Cha 6 


2003 Jan 20 


S27 


6x16 


6x11 


6x11 


VIS-VIS 


1.3 


110 


93 


75 


1 .65 




j 5 


22 


Sz 48 


2004 Apr 06 


S13 


4x48 






IR-K 


0.6 


184 






1 .66 








BK Cha 


2004 Apr 03 


S13 


4x30 






IR-K 


0.5 


78 






1 .66 


-j 






Sz 60 


2004 Apr 05 


S13 


4x60 






IR-K 


0.5 


115 






1 .66 








Sz 62 


2003 Jan 22 


S27 


15x3.0 


15x1.6 


15x1.6 


VIS-VIS 


1.0 


150 


105 


91 


1.68 


3 


7 


11 


Herschel 4636 .... 


2003 Jan 20 


S13 


9x3.0 


9x4.5 


9x4.5 


VIS-VIS 


1.0 


70 


1 12 


104 


1.08 


1 s 


] 7 


17 


ESO Hff 283 


2003 Mar 26 


S13 


3x46 


3x35 


3x33 


VIS-VIS 


0.7 


76 


74 


75 


1 .36 


13 


35 


40 


Sz 65 


2004 Apr 06 


S13 


36x2.2 






VIS-VIS 


0.8 


60 






1 .36 


20 






Sz 68 


2004 Apr 06 


S13 


88x0.6 






VIS-VIS 


0.6 


61 






1.34 


27 






HO Lup 


2004 Apr 06 


S13 


4x13.8 






IR-K 


0.7 


80 






1.37 


7 






Sz 101 


2004 Apr 10 


S13 


4x18 






IR-K 


0.7 


79 






1 .03 


10 






Sz 108 


2004 Apr 10 


S13 


8x8.0 






IR-K 


0.6 


76 






1 .03 


13 






Sz 120 


1AA.1 A T A 

2004 Api 1 


S13 


180x0.6 






WTO WTO 

VIS-VIS 


0.7 


70 






1 .04 


18 






WSB 3 


2004 May 01 


S13 


4x16.8 






IR-K 


0.9 


133 






1.10 


4 






WSB 11 


2004 Jun 18 


S13 


4x56 






IR-K 


1.1 


93 






1.04 


2 






WSB 20 


2004 May 01 


S13 


12x5.0 






IR-K 


0.8 


75 






1 .06 


12 






WSB 28 


2004 May 01 


S13 


4x24 






IR-K 


0.8 


117 






1.04 


8 






SR 24 


2004 May 01 


S13 


16x4.0 






IR-K 


1.1 


72 






1 .03 


20 






Elias 2-30 


2004 May 01 


S13 


40x2.0 






IR-K 


1.0 


72 






1.02 


13 






WSB 46 


2004 May 01 


S27 


4x13.6 






IR-K 


0.8 


148 






1.01 


3 






Haro l-14c 


2004 May 01 


S27 


52x1.0 






VIS-VIS 


1.1 


82 






1.07 


7 






ROXs 43 


2004 May 04 


S13 


100x0.6 






VIS-VIS 


0.5 


80 






1 .06 


36 






Elias 2-49 


2004 May 01 


S13 


56x1.2 






VIS-VIS 


0.9 


77 






1 .03 


1 1 






HBC 652 


2004 May 01 


S27 


100x0.5 






IR-K 


0.9 


84 






1 .03 


9 






LkHff 346 


2004 May 01 


S27 


16x4.0 






IR-K 


0.9 


104 






1.02 


16 







a : The seeing value quoted is a median value of all frames at 550 nm at zenith. 
h : Observed with a neutral density filter. 

c : Only 2 acquisition frames in Bry with very poor Strehl ratio. 



3. Observations and data reduction 

Observations were carried out during two periods. A first set of 
37 objects were observed from October 22' /! 2002 to March 26'' ! 
2003 while observations of another 21 systems were conducted 



from April 4" ! 2004 to June \l' h 2004. Magnitudes from the lit- 
erature in V and from 2MASS (Two Micron All-Sky Survey) 
in K are reported in Table^for each combined system or alter- 
natively for the brightest components if resolved. NACO was 
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5 10 15 

separation [arcsec] 

Fig. 1. Distribution of separations/orbital semimajor axes of the 
wide binaries of our sample. 



mainly used with the visible wavefront sensor (WFS) since 
most of the objects resolved by the wavefront camera and used 
as natural guide star for the AO system have bright enough V- 
magnitude. Only in a few cases the use of the IR-WFS was 
needed. Each object of the first set has been observed through 
the three narrow-band filters Bry (2.166yum, 0.023/mi width), 
H 2 (2. 1 22 /an, 0.022 fim width) and [Fell] ( 1 .644 jjm, 0.0 1 8 /im 
width) in order to prevent saturation from the brightest stars 
of our sample. Objects of the second set were observed only 
through the [Fell] filter. Three or four dithered exposures were 
taken in each filter. Single exposure integration times vary from 
typically 0.5 s to 30 s which imply a magnitude limit for the 
median exposures of H ~ 14 and K ~ 13, and for the deepest ex- 
posures H ~ 16 and K ~ 15.5, respectively (see Tables[5]and0}. 
The field of view (FOV) is either ~ 13'.'6xl3'.'6 or ~ 27'.'8x27'.'8 
for the S13 and S27 cameras respectively, with some enlarge- 
ment which can be obtained with the dithering depending on 
the position of the target in each frame. The combination of 
natural guide star magnitude and seeing lead to AO-corrections 
with typical Strehl ratios of ~30% and ~ 10% in Bry and 
[Fell], respectively, which provides mainly diffraction-limited 
cores. The resolutions of the images, measured as the FWHM 
of a Gaussian fit to the radial profile of the brightest compo- 
nent, have a median value of 80.2 mas, 79.7 mas and 77.8 mas 
in [Fell], H2 and Bry, respectively (see Tables- 
Data reduction was performed in IRAF 1 in the usual way : 
sky subtraction, flat-fielding, bad-pixels and cosmics correc- 
tions. Sky background frames obtained from median averag- 
ing of the dithered frames were subtracted from the individuals 
frames. All frames were then flat-fielded and corrected for bad- 
pixels using calibration files provided by ESO. Cosmic rays 
have been removed using the LACOS 2 package (Van Dokkum 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract to the National Science 
Foundation. 

2 LACOS webpage : http://www.astro.yale.edu/dokkum/lacosmic/ 



2001 ) and finally the frames were registered and combined to- 
gether. 

The detection of companions was first performed by eye. 
We then searched for eventual close and faint companions 
using psf-subtraction. Specifically, each component image of 
each system was subtracted by the other scaled component im- 
ages of that system and the residuals analyzed. The efficiency 
of this method is substantially improved by the fact that the 
components were imaged simultaneously and that their sep- 
aration are usually within the isoplanatic angle, i.e. the AO- 
correction is very similar. No companion undetected by eye 
was found, however a few components exhibit residuals which 
are presumably to be attributed to faint extended emission of 
scattered light from disks or envelopes (an example is presented 
in AppendixlAl. 

Astrometry and relative photometry of the high-order 
multiple systems were performed using the IRAF package 
DAOPHOT. The relative photometry was based on the results 
of PSF-fitting using the task ALLSTARS (Stetson [19871 aver- 
aged with the results obtained with aperture photometry tech- 
niques, with radius of typically 3 pixels. The PSF used was 
usually the main component of each system. Relative positions 
were determined with the PSF-fitting. To estimate uncertain- 
ties, we performed the same analysis on both the individual 
frames and the combined frames for all wavelengths and cal- 
culate their standard deviation. We assumed that the system- 
atics introduced by the PSF-fitting procedure are of the or- 
der of 07001 and 0.01 in the derived separation and flux ra- 
tio, respectively. The final error in relative positions are es- 
timated by combining quadratically the rms variations of the 
results from the different data sets with the uncertainty in the 



plate scale (S13 : 13.26+0.03 mas, Masciadri et al. l2003l S27 : 
27.01 +0.05 mas, Chauvin et al. 2005) and detector orientation 



(±0.5°, Masciadri et al. 2003 ). The astrometric and photometric 
results are presented in Table[3]for all possible triple/quadruple 
systems of our sample. 

4. Results 

4.1. New candidate companions 

As a result of our survey, among the 58 observed systems, 
10 were found to be apparent triples and 7 to be apparent 
quadruples (see Fig. [2] and Table[3j. The remaining systems 
remain binary within our resolution and detection limits (see 
Table0}. No new edge-on disks other than the two previously 
known in our sample (HK Tau B and LkHa 263 C) were de- 
tected. However, since LkHa263 C was not detected it is prob- 
able that our survey might not have been deep enough to detect 
those (see AppendixlBl. Throughout this paper, the single com- 
ponents of all systems are named with upper case letters, start- 
ing from A for the brightest component in Bry and [Fell], for 
respectively the first and second data set. This nomenclature 3 is 
adopted in Fig. [2 Our designation is only practical when deal- 
ing with multiple systems, and does not reflect primary, sec- 
ondary, tertiary etc. which are based on the component masses. 



3 It should be noted that the Multiple Star Catalogue (MSC, 
Tokovinin 1997 1 uses a more sophisticated nomenclature. 
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CGHa5 



COH.S/6 
Bry 



CG Ha 6 

A B 



ROXs 43 
Fell 



Sz 68 
Fell 



l.kH. 262/263 
Bry 



SR24 
Fell 



UZTan 
Bry 



Ud 1.336 
Bry 



Sz 30 
Bry 



PH„ 30 

• Bry 

A 



1 1 1 I 1 — 



Sz.4l 
Bry 



Lkll.346 
Fell 



J 4872 
Bry 



* B 
♦ 

. C 



KSO H„ 283 
Bry 



DX Tau 
Bry 



A 


• 


B 





HV 1 an 
Bry 



VWCha 
Bry 



B GO Tau 

* Bry 



Fig. 2. Already known and candidate triple (upper panel) and quadruple (lower panel) systems detected in our VLT/NACO 
survey. These systems are the only apparent high-order multiple systems found in our sample, physical systems will be identified 
in Sect. 14.51 North is up, east is left. 
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Table 3. Measured parameters of all candidate triple/quadruple systems resolved in our survey, with their apparent multiplicity 
in the projected separation range 0707-12". Also quoted in last column is the probable and assumed multiplicity of this systems 
in this separation range based on the chance projection analysis of Sect. 14.51 



System 


app. 


Pair 


Separation 


P.A. 




Flux ratio 




1 st component pair 


2nd component pair 


Ref 


Prob. 




mult. 




[arcsec] 


[°] 


[Fell] 


H 2 


Bry 


H 


K 


H 


K 


(mult.) 


mult. 


ESO Ha 283.... 


Q 


A-B 

A-C 
A-D 


2.058 ±0.002 
2.571 ±0.003 
y.jzo ± u.uu / 


238.5 ±0.5 
165.1 ±0.5 

254.6 ± 0.5 


0.61 ±0.01 
0.04 ±0.01 

\)X)5 ± U.U1 


0.42 ±0.01 
0.06 ±0.01 
U.UJ ± U.Ui 


0.42 ±0.01 
0.06 ±0.01 
U.Uj ± w.wl 


11.64 ±0.06 


10.69 ±0.03 


12. 18 ±0.07 
15.10±0.17 
14.77 ±0.15 


11.64 ±0.05 
13.79±0.07 
14.09 ±0.10 


1 


B 


I'll.) 30 


T 


A-B 

A-C 


0.661 ±0.002 
11.58±0.01 


350.3 ±0.5 
195.7 ±0.5 


0.25 ± 0.01 
0.99 ±0.02 


0.25 ±0.01 
1.02 ±0.02 


0.24 ±0.01 
0.95 ± 0.01 


12.70 ±0.03 


12.39 ±0.03 


14.19 ±0.06 
12.81 ±0.02 


13.97 ±0.05 
12.60 ±0.03 


1 


B 


CGHo- 5/6 


T 


A-B 
A-C 


0.482 ±0.001 
11.289 ±0.007 


306.8 ±0.5 
317.4±0.5 


0.59 ±0.01 
0.31 ±0.03 


0.91 ±0.01 
0.54 ±0.03 


0.92 ±0.01 
0.60 ±0.01 


10.04 ±0.02 


9.82 ±0.03 


10.61 ±0.03 
10.82 ±0.03 


9.91 ±0.03 
10.02 ±0.02 


1 


T 


VW Cha 


Q 


A-B 
B-C 
A-D 


0.653 ±0.001 
0.113±0.001 
16.78 ±0.01 


177.2 ±0.5 
232.9 ±0.5 
220.4 ±0.5 


0.39 ±0.01 
0.62 ±0.05 
0.03 ±0.01 


0.25 ±0.01 
0.61 ±0.03 
0.04 ±0.01 


0.24 ±0.01 
0.62 ±0.04 
0.05 ± 0.01 


8. 17 ±0.07 


7.32 ±0.04 


9.19±0.10 
9.71 ±0.19 
10.86 ±0.02 


8.87 ±0.07 3, 9, 10 
9.40±0.13 
9.89 ±0.02 


Q 


Sz 30 


T 


A-B 
B-C 
A-C 


1.237 ±0.004 
3.854 ±0.004 
4.504 ±0.002 


19.4 ±0.5 
85.3 ±0.5 
70.8 ±0.5 






0.77 ±0.01 
0.86 ±0.05 
0.66 ±0.03 




9.58 ±0.05 




9.86 ±0.05 
9.98 ±0.03 


2 


T 


Sz41 


T 


A-B 
A-C 


I. 974 ±0.002 

II. 46±0.01 


162.4 ±0.5 
66.1 ±0.5 


0.13±0.01 
0.23 ± 0.02 


0.11 ±0.01 
0.26 ±0.02 


0.12 ±0.01 
0.22 ±0.02 


8.66 ±0.06 


8. 12 ±0.03 


10 87 ± 07 
8.72 ±0.05 


10.45 ±0.04 
8.57 ±0.02 


2, 15 


B" 


LkHti- 336 


T 


A-B 
B-C 
A-C 


10.777 ±0.008 
8. 199 ±0.007 
5.799 ±0.004 


95.5 ±0.5 
243.4 ±0.5 
144.2 ±0.5 


0.26 ±0.01 
0.46 ±0.01 
0.12±0.01 


0.27 ±0.01 
0.39 ±0.01 
0.10±0.01 


0.27 ±0.01 
0.41 ±0.01 
0.11 ±0.01 


9.72 ±0.03 


9.15 ±0.03 


12.00 ±0.05 
11.41 ±0.03 


11.47 ±0.03 
10.73 ±0.02 


2 


T 


J 4872 


Q 


A-B 
CD 

A-C 


0.173 ±0.001 
0.102 ±0.001 
3.339 ±0.003 


74.8 ±0.5 
157.9 ± 1.3 
234.1 ±0.5 


0.86 ±0.01 
0.76 ± 0.04 
0.44 ± 0.04 


0.86 ±0.01 
0.71 ±0.01 
0.49 ±0.01 


0.85 ± 0.01 
0.73 ±0.01 
0.50 ±0.01 


9 53 ±0 08 
10.42±0.17 


9.63 ± 0.05 
10.39 ±0.06 


9.69 ±0.09 
10.71 ±0.14 


9.80 ±0.05 
10.74 ±0.07 


4,5, 1 


Q 


UX Tau 


Q 


A-B 
B-C 
A-D 


5.856 ±0.003 
0.1 36 ±0.001 
2.692 ±0.002 


269.7 ±0.5 
309.0 ± 0.5 
181.6±0.5 


0.12±0.01 
0.80 ±0.01 
0.07 ±0.01 


0.11 ±0.01 
0.81 ±0.02 
0.05 ±0.01 


0.13 ±0.01 
0.81 ±0.02 
0.05 ±0.01 


8.03 ±0.02 


7.61 ± 0.02 


9.58 ± 0.05 
9.83 ± 0.03 
10.94 ±0.03 


9.56 ±0.02 
9.80 ±0.04 
10.79 ±0.04 


2,4 


Q 


UZ Tau 


T 


A-B 
B-C 


3 560 ± 006 
0.365 ±0.006 


273 5 ±0 5 
4.3 ±3.0 


22 ± 04 
0.52 ±0.03 


22 ± 02 
0.51 ±0.01 


24 ±0 01 
0.52 ±0.01 


8.12 


7.35 ± 0.03 


8.46 ±0.04 
9.18 ± 0.11 


7.93 ±0.03 
8.64 ±0.05 


6, 7, 8 


T 


FV Tau 


Q 


A-B 
CD 
A-C 


0.780 ±0.004 
0.693 ±0.003 
12.081 ±0.009 


273.2 ±0.6 
294.6 ±0.7 
106.1 ±0.5 


0.16 ±0.01 
0.11 ±0.02 
0.71 ±0.01 


0.24 ±0.01 
0.15 ±0.02 
0.56 ±0.01 


0.29 ±0.01 
0.18 ±0.02 
0.51 ±0.01 


8.49 ±0.04 
9.60 ±0.04 


7.72 ± 0.03 
9.05 ± 0.04 


10.48 ±0.11 
12.03 ± 0.20 


9.07 ±0.08 
10.91 ±0.13 


6 


Q 


GG Tau 


Q 


A-B 
CD 

A-C 


0.249 ±0.001 
1.460 ±0.003 
10. 100 ±0.007 


346.0 ±0.5 
135.3 ±0.5 
186.0 ±0.5 


0.58 ± 0.02 
0.21 ±0.01 
0.08 ± 0.01 


0.57 ±0.01 
0.20 ±0.01 
0.08 ±0.01 


0.58 ±0.01 
0.21 ±0.01 
0.08 ± 0.01 


8.31 ±0.04 
10.60 ±0.03 


7.86 ±0.02 
10.18 ±0.03 


8.91 ±0.07 
12.30 ±0.06 


8.46 ±0.04 
11.87 ±0.07 
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Q 


LkHo- 262/263... 


Q" 


A-B 
A-D 


0.414 ±0.001 
15.38 ±0.01 


52.0 ±0.5 
202.5 ± 0.5 


1.12 ±0.02 
0.90 ±0.06 


0.87 ±0.01 
1.05 ±0.03 


0.88 ±0.01 
0.99 ±0.05 


10.66 ±0.03 


10.18 ±0.02 


10.53 ±0.04 
10.29 ±0.02 


10.32 ±0.03 
9.64 ±0.02 


12, 13 


Q 


LkHtr 346 


Q 


A-B 
A-C 
A-D 


0.204 ±0.002 
1.086 ±0.001 
5.068 ±0.006 


256.5 ± 1.0 
328.2 ±0.5 
329.7 ±0.5 


0.24 ±0.03 
0.26 ±0.01 
0.60 ±0.01 






9.23 ±0.07 




10.8 ±0.2 
10.70 ±0.08 
8.99 ±0.03 




1 


Q 


ROXs 43 


T 


A-B 
A-C 


0.334 ±0.001 
4.487 ±0.001 


168.9 ±0.5 
12.1 ±0.5 


0.05 ± 0.01 
0.64 ±0.01 






7.22 ±0.06 




10.54 ± 0.07 
7.48 ± 0.08 




1 


T 


SR 24 


T 


A-B 
B-C 


5.068 ±0.001 
0.081 ±0.001 


348.5 ±0.5 
45.6 ±0.5 


0.17 ±0.01 
0.5 ±0.1 






8. 17 ±0.04 




9.1 ±0.2 
9.8 ±0.4 




14 


T 


Sz 68 


T 


A-B 
A-C 


0.126 ±0.001 
2.808 ±0.002 


246.7 ±0.5 

296.8 ±0.5 


0.16±0.01 
0.07 ±0.01 






7.09 ±0.03 




9.05 ± 0.05 
10.05 ±0.05 




9 


T 



a : includes the edge-on disk LkHo' 262/263C (Jayawardhana et al. 2002 Chauvin et al. 2002 i undetected in our survey. 
b : Sz41 C was unambiguously identified as a background giant fWalte ri!992l . 
Notes : multiplicity abbreviations are B=Binary, T=Triple and Q=Quadruple. 

Multiplicity r eference s : (1) This work. (2) RZ93 . (3) Brandeker et alJOTfl (4) Duchen dTT3991 (5 ) White & GhezJSM] (6) Simon et al fiWI (7) L einert e t al fiWfl ( 8) Ghe z et al JlWfl 
(9) Ghez et al. 1997a (10) Brandner et al. 1996" (1 1) Leinert et al J199ll (12) Jayawardhana et al j2TOI (13) Chauvin et al j2TOI (14) Simon et al JT995l (15) Walter [TWI 



We mainly searched for unresolved pairs in each individual 
component of the wide binaries of the sample. Nevertheless 
some of the targets of the sample were so poorly known and 
observed that it has been possible to detect new distant com- 
panion candidates. Indeed there is, in some cases, an obvious 
lack of NIR follow-up observations since the survey of RZ93 
which was performed with Gunn-z filter (A ~Q.9fim). 

Three of the triple candidate (CGHa5/6, PHa30 and 
ROXs 43) and two of the quadruple candidate systems (ESO 
Ha 283 and LkHo- 346) are new detections. CGHa6 was 
known from RZ93 to form a binary with the 11'.' 1 north- 
western companion CGHa5. Here we resolved for the first 
time CGHa6 as a close 075 binary, making it a hierarchical 
triple. PHa30 is a 0767 binary with a wide companion found 
in our survey at 11'.'9 separation. ESO Ho- 283 appears to be a 
possible quadruple system composed of the original 2'.'08 bi- 
nary detected by RZ93 plus a third (C) and a faint fourth (D) 
component located 2'.'57 SE and 973 SW of the primary, re- 



spectively. The LkHo- 346 system contains two additional can- 
didate companions. The first one is located ~072 from the pri- 
mary of the system (spectral type K7, aka LkHa346NW), 
while the second one is lying approximately in between the pri- 
mary and the secondary (spectral type M5, aka LkHa-346 SE). 
With flux ratios of ~0.3 with respect to the M5 component, 
both might well be substellars, if bound. This system is further 
discussed in AppendixE] ROXs 43 A (NTTS 162819-2423S) 
is a 89-days single-lined spectroscopic PMS binary in the p 
Oph cloud (Mathieu et al. 1989 1, with eccentricity e=0.4 and 
a.sin(i)=0.1 AU. It is part of a hierarchical quadruple system, 
ROXs 43 A lying 4'.'8 south from R OXs 43 B (NTTS 162819- 
2423N), a 07016 binary (Simon et al. [T995l . We found a candi- 
date fifth component in this system 073 from ROXs 43 A, corre- 
sponding to 42 AU at 140 pc, with flux ratio ~0.05 in H-band. 
Interestingly, this candidate companion was not detected in a 
recent speckle survey of p Oph (Ratzka et al. 2005 ). Part of the 
near-IR excess of ROXs 43 A could originate from this candi- 
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Table 4. Upper limits for the relative brightness of undetected companions to both components of the binaries in our survey, 
measured at 0707, 0715 and 075, as well as the limiting magnitudes at large separations. 



Upper limits in [Fell] Upper limits in Bry Limiting 



System 


Separation 


PA." 


Flux ratio a 




comp. A 




comp. B 




comp. A 




comp. 


B 


Magnitudes 




[arcsec] 


[°] 


[Fell] 


Bi y 


0707 


O. 15 


0. 


0707 


u. 15 


0. ^ 


0707 


0. 15 


0. 5 


0707 


(X' 15 


075 


H 


K 


AR Ori 


2.002 ±0.001 


239.4 ±0.5 


74 


55 


0.19 


18 


08 


0.18 


17 


o i 


0.17 


12 


05 


0.17 


0.12 


0.1 1 


1 3 S 


1 3 6 


CoD -29°8887 . 


0.515 ±0.003 


28.2 ±0.5 


0.52 


5 1 


0.04 


(p 


0.007 


0.04 


0.02 


0.02 


0.03 


0.01 


0.005 


0.03 


0.01 


0.01 


P 9 


|9 x 


CO Ori 


^ i\CO i A I\t\A 

2.05ft ± 0.004 


274.7 ± 0.5 


0.08 


0.07 


0.39 


0.13 


0.02 


0.39 


0.31 


25 


0. 15 


0.1 1 


0.007 


0.15 


0.14 


0.11 


12.3 


1 1 * 


CV Cha 


i i oai _i_ a nn< 
1 1 .ZOj ± O.OOj 


98.7 ± 0.5 


1 1 


14 


0.76 


2^ 


04 


0.79 


1 1 6 i 


5 1 


0.26 


14 


02 


0.28 


0.30 


0.23 


111 


1 1 7 


DK Tau 


Z.joO ± 0.001 


1 1 A A _i_ A ^ 
1 lv.0 ± 0.5 


i ^ 


7 s 


0.06 


05 


1 


0.06 


04 


03 


0.04 


04 


1 


0.05 


0.04 


0.03 


1 3 4 


1 1 9 


eju tia zoi 


1 1 AO i A A A 1 
1 . /4S ± 0.001 


172.6 ± 0.5 


1 1 


7 2 


0.09 


0.07 


0.02 


0.09 


0.07 


0.06 


0.09 


0.05 


()| 


0.09 


0.05 


0.05 


1 5.3 


1 5. 1 


Glass l b 


2.430 ± 0.002 


285.1 ±0.5 


0.27 


0.23 


0.10 


0.07 


0.02 


0.09 


0.06 


0.05 


0.05 


0.02 


0.004 


0.04 


0.02 


0.02 


13.0 


1 3.4 


HD 76534 


2.073 ±0.002 


303.6 ±0.5 


0.36 


0.33 


0.09 


0.08 


0.02 


0.09 


o.os 


0.07 


0.07 


0.07 


0.02 


0.07 


0.07 


0.07 


12.4 


1 2.4 


Hen 3-600 


1.481 ±0.003 


213.4 ±0.5 


0.69 


0.68 


0.05 


0.04 


0.02 


0.05 


0.04 


0.03 


0.03 


0.03 


0.01 


0.03 


0.03 


0.02 


1 1 .8 


1 2.2 


Herschel 4636 ... 


3.754 ±0.004 


213.8 ±0.5 


0.23 


0. 1 1 


0.10 


0.05 


0.005 


0.10 


0.05 


0.03 


0.42 


0.08 


0.004 


0.42 


0.09 


0.02 


14.0 


1 3.6 


HK Tau 


2.318 ±0.006 


171.3 ±0.5 






0.70 


u.uo 


U.U.I 








0.61 


U.U4 


u.uz 








1 J.Z 


I J.Z 


HN Tau 


3.142±0.001 


219.7 ±0.5 




0.12 














0.31 


0.09 


0.02 


0.41 


0.33 


0.33 




I 3.0 


IT Tau 


2.416 ±0.008 


225.1 ±0.5 


28 


7 7 


0.40 


2? 


07 


0.39 


36 


30 


0.14 


1 2 


(I (p 


0.15 


0.14 


0.13 


12 


1 7 y 


L1642-1 


O 1AA _i_ A AA1 
Z. /44 ± 0.001 


j40. / ± O.J 


36 


7 3 


0.17 


n l s 

U. 1 J 


U.UJ 


0. 17 


15 


| ^ 


0. 1 1 


n no 
u.uv 


u.uz 


0. 1 1 


0.09 


0.09 


12 5 


1 "* 


LKriO' Z06 


1 T)/l _i_ A (\(Yl 

J.ZZ4 ± O.OOj 


1 in T j A c 
1 / /. / ± O.J 


0.79 


0.95 


0.10 


0.05 


0.02 


0. 1 


() ()S 


0.03 


0.07 


0.04 


( p 


0.07 


0.04 


0.03 


14.0 


13.9 


PHa 14 


A £CA i A AA£ 

0.650 ± 0.006 


168.3 ± 0.5 


85 


9 i 


0.05 


05 


DO 7 


0.05 


04 


o 7 


0.06 


05 


(i (p 


0.06 


0.05 


0.02 


16 1 


15 7 


PHtv 51 


O.oZ / ± O.OOj 


1 OA A i A C 

184.0 ± O.J 


59 


46 


0. 19 


18 


07 


0. 1 9 


18 


| ^ 


0. 17 


16 


05 


0. 1 8 


0. 16 


0. 13 


1 5 7 


14 ^ 


mi; a 

KW Aur 


i a a n i a nni 
1.44/ ±0.003 


255.0 ± 0.5 


0^20 


0^21 


0.14 


0'08 


o'm 


0. 1 6 


0.14 


o'.w 


0.08 


(106 


0.007 


0.08 


0.06 


0.04 


12.9 


13^1 


SX Cha 


Z.Z1 1 ± O.OOj 


TOO 1 i A C 

Z80. / ± O.J 


0.25 


0.22 


0. 12 


0.09 


0.01 


0. 1 3 


0.09 


0.05 


0. 13 


0.11 


0.02 


0. 13 


0. 1 1 


0. 10 


15.3 


13.5 


Sz 15 


10.506 ±0.008 


27.5 ±0.5 


0.45 


0.60 


0.49 


0.20 


0.03 


0.53 


0.31 


0.12 


0.44 


0.14 


0.02 


0.50 


0.29 


0.09 


14.9 


14.9 


Sz 19 










0.43 


0.07 


0.04 








0.56 


0.07 


0.04 








10.5 


9.7 


Sz 62 


1.1 16 ±0.005 


267.3 ±0.5 


o!eo 


0'59 


0.09 


0.09 


0.03 


oii'o 


0i)8 


o!o6 


0.08 


0.07 


0.03 


01)8 


om 


0.06 


14.2 


13.7 


vBH 16 


4.509 ±0.003 


329.8 ±0.5 


0.59 


0.60 


0.20 


0.18 


0.09 


0.21 


0.18 


0.13 


0.14 


0.10 


0.04 


0.14 


0.10 


0.06 


12.5 


12.8 


VV Cha 


0.790 ±0.001 


13.4 ±0.5 


0.72 


0.63 


0.22 


0.19 


0.10 


0.21 


0.19 


0.14 


0.25 


0.23 


0.10 


0.24 


0.23 


0.17 


13.1 


12.6 


BK Cha'' 


0.766 ±0.001 


153.6±0.5 


0.55 




0.04 


0.02 


0.01 


0.05 


0.02 


0.01 














15.4 




Elias 2-30 


6.317 ±0.001 


175.3 ±0.5 


0.06 




0.28 


0.1 1 


0.005 


0.28 


0.13 


0.10 














13.7 




Elias 2-49 


1.097 ±0.001 


224.9 ±0.5 


0.17 




0.04 


0.02 


0.005 


0.04 


0.02 


0.01 














13.7 




Haro l-14c 


13.047 ±0.001 


122.6 ±0.5 


0.42 




0.28 


0.03 


0.003 


0.92 


0.27 


0.06 














14.4 




HBC 652 


8.695 ±0.001 


321.7±0.5 


0.13 




0.23 


0.07 


0.005 


0.18 


0.07 


0.02 














14.7 




HO Lup 


1.536±0.001 


34.9 ±0.5 


0.24 




0.09 


0.06 


0.02 


0.09 


0.06 


0.05 














14.3 




Sz 101 


0.771 ±0.001 


302.9 ±0.5 


0.54 




0.04 


0.02 


0.01 


0.05 


0.02 


0.02 














15.3 




Sz 108 


4.029 ±0.001 


24.5 ±0.5 


0.14 




0.07 


0.04 


0.009 


0.07 


0.04 


0.04 














14.8 




Sz 120 


2.600 ±0.001 


140.8 ±0.5 


0.09 




0.03 


0.01 


0.002 


0.03 


0.01 


0.01 














14.1 




Sz 48 


1.480 ±0.001 


232.5 ±0.5 


0.58 




0.12 


0.10 


0.06 


0.12 


0.10 


0.08 














14.2 




Sz 60 


3.347 ±0.001 


99.7 ±0.5 


0.80 




0.07 


0.05 


0.02 


0.07 


0.05 


0.02 














14.6 




Sz 65 


6.371 ±0.001 


97.8 ±0.5 


0.31 




0.06 


0.03 


0.005 


0.06 


0.03 


0.02 














14.5 




WSB 11 


0.472 ±0.001 


316.8 ±0.5 


0.89 




0.11 


0.09 


0.04 


0.11 


0.09 


0.06 














15.2 




WSB 20 


0.779 ± 0.002 


31.6±0.5 


0.46 




0.02 


0.01 


0.006 


0.02 


0.01 


0.01 














14.5 




WSB 28 


5.119±0.001 


358.3 ±0.5 


0.07 




0.27 


0.26 


0.03 


0.27 


0.27 


0.23 














14.0 




WSB 3 


0.600 ±0.004 


160.2 ±0.5 


0.65 




0.08 


0.07 


0.03 


0.08 


0.07 


0.06 














14.0 




WSB 46 


10.209 ±0.001 


301.6 ±0.5 


0.50 




0.24 


0.16 


0.03 


0.24 


0.15 


0.04 














14.6 





a : Assumes as for the triples/quadruples that the primary is the brightest component in Bry and [Fell], for respectively the first and second data set. 
b : Components inverted with respect to RZ93, given the above assumption. 

Notes : Separations and PA.s are computed from aperture photometry and the uncertainties combine the rms variations over the individual frames and the uncertainties due to the plate 
scale and detector orientation. 



date companion. Jensen & Mathieu ( 1997 > fitted a circumbi- 
nary disk with gap to the SED of ROXs 43A. If the candidate 
companion is physically related, the issue of the spatial distri- 
bution of circumstellar/cicumbinary material would have to be 
re-addressed. Additionally, J 4872 B was already known to be 
a close pair (White & Ghez-2001f, forming a quadruple with 
the other pair J 4872 A, but no measurements were published 
before ours. 

4.2. Limits for undetected companions and 
completeness 

Table|4] shows the upper limits for undetected companions to 
the only binaries of our sample at separations 0707, 0715 and 
075 in [Fell] and Bry, as well as the estimated limiting mag- 
nitudes at large separations in H and K. The detection limits 
were computed using the following method. At each radial dis- 
tance and position angle from a star, the standard deviation of 
the flux was calculated over a circular region of radius 70 mas, 
i.e. equivalent to the mean size of the PSF cores. The detection 
limit as a function of separation to a star is the average of the 
5<t flux over all position angles excepted those lying in the di- 
rection of the companion. We repeated the procedure with the 



PSF-subtracted images, and defined the total detection limit as 
the minimum of the two detection limits. The sensitivity is in- 
creased by the PSF-subtraction at separations smaller than typ- 
ically 071-073. Our procedure was checked by adding artificial 
companions at a few separations and position angles, which 
are scaled versions of the target star PSF, and looking for the 
maximum flux ratio at which the companion is detected by vi- 
sual inspection. A statisfying agreement was found between the 
limiting flux ratio given by the two methods. 

Concerning the limiting magnitude at large separations, 
these were calculated by combining 2MASS H- and K- 
magnitudes of resolved or unresolved systems and our limit- 
ing flux ratios in [Fell] and Bry. We checked our estimated 
magnitude using undetected or marginally detected stars in our 
frames which are in common with 2MASS and found a rel- 
atively good agreement. The discrepancies are up to ~lmag 
and can be explained either by possible photometric variabil- 
ity between 2MASS and our observation epochs or by the fact 
that we assumed that the relative flux in H and K is given by 
that in [Fell] and Bry. The more remarkable cases are further 
discussed in Sect. 14. 51 

For any given component of the binaries we thus have the 
limiting flux ratio for undetected companion as a function of 
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Fig. 3. Relative motion in RA and Dec between the RZ93 po- 
sitions and those measured in this study (NACO-RZ93). Green 
triangles are pairs with projected separations < 150 AU. 

separation. We can therefore produce a completeness map for 
each of these sources, i.e. a map giving the probability to de- 
tect a companion as a function of separation and magnitude dif- 
ference for each source. The underlying assumption is that all 
companions above the sensitivity limits are detected, i.e. there 
is a sharp limit in sensitivity. This might not be completely true 
especially for AO imaging data because of the anisotropy of 
the correction, but this simplification should still provide mean- 
ingful results (see Tokovinin et al. 2006 and Kohler et al. 2006 
for further discussions). The total completeness map of all the 
sources is simply the average of all individual completeness 
maps. We will further assume in the following that the total 
completeness map created this way is statistically equivalent to 
that obtained if we would have included the triple and quadru- 
ple systems. Fig.@]shows the total completeness map in H and 
K with the detected companions of the triple/quadruple systems 
overploted. The majority of the companions fall above the 90% 
completeness level and all are above the 50% level, giving us 
confidence in these limits. 

Assuming that both the distribution of separation and the 
distribution of flux ratio of companions are flat, which is obvi- 
ously a very rough assumption, we estimate that the complete- 
ness within the separation range 0707-10" and above a flux 
ratio of 0.001 is of the order of 95%. Given the statistical er- 
rors involved in this study, which are a direct consequence of 
the sample size, the correction for incompleteness is considered 
not significant enough to be applied. It should be noted that sys- 
tems observed with the S13 camera were not probed for com- 
panions at large separations (says 5; 7-8", depending of the po- 
sition of the system on the detector). Even though the brighter 
companions at such separations would show up in 2MASS, the 
faintest would likely be missed. 

4.3. Microjets and Herbig-Haro objects 

Since the exposure times for all sources was scaled in order to 
avoid saturation, the sensitivity to extended emission is lim- 



ited. We could however be sensitive to bright Herbig-Haro 
(HH) objects and/or HH-jets close to the star by using PSF- 
subtraction techniques in H2 and [Fell]. This is of particu- 
lar interest in view of the known correlation of outflows and 
multiple stars (Reipurth 2000), even associated with relatively 
old sources like classical T Tauri stars (e.g. Mundt & Eisloffel 
1998). Similarly, dynamical interaction in multiple stars sys- 
tems following an ejection of one component by close triple en- 
counter are also suspected to be at the origin of the FU Orionis 
phenomenon (Reipurth & Aspin 2004). 

Our clear non detection of any prominent outflows/jets 
could however be explained by a lack of sensitivity since some 
sources known to harbor jets were not detected. The bipolar 
jet originating from RW Aur A (HH 229) is a good example 
in this respect. First identified spectroscopically by Hirth et al. 
( 1 994 ), its inner part (07 1-10") has later been imaged as a bipo- 
lar micro-jet in the optical forbidden emission-lines [OI] and 
[SII] (Dougados et al. 2000 1, but was already known to have 
a large extension, up to 100" (Mundt & Eisloffel 1998), and 
could be even more extended if HH 835, some 5' distant, is as- 
sociated with the outflow (McGroarty & Ray 2004). Recently, 
the microjet has been detected in [Fell] (Davis et al. 2002 Pyo 
et al. 2005 ), and both high and low-velocity components iden- 
tified. Our detection limit is H=12.9 which translates to ~9.4 
mag.arcsec -2 , which is probably l-2mag brighter than typical 
[Fell] fluxes measured from HH-objects (e.g. Gredel I1994> . 
Similarly, probably for the same reason, the recently discov- 
ered microjet in UZ Tau E, imaged in [OI] using HST/STIS 
slitless spectroscopy (Hartigan et al. 2004 1, was not detected in 
our images in [Fell] and H2. 

However, we note that we might have detected some out- 
flow emission from another known HH-jet, namely HH 186, 
whose driving source is Sz 68 (Heyer & Graham 1989). This is 
a highly collimated jet of forbidden [SII] emission extending 
34" at a position angle of 135°, and composed of three emis- 
sion peaks at 21", 28" and 34". These knots are outside our 
field of view and we did not detect any jet emission on this 
side of Sz 68 A. But, we detected a faint nebulosity at PA~295° 
and separation ~ 1727 which could be the possible counter jet 
quoted in Heyer & Graham ( 1989). Deeper imaging is neces- 
sary in order to confirm this result. 

4.4. Relative motions 

Fig-El shows the variations of relative positions between our 
measurements and the values reported in RZ93. Uncertainties 
of the RZ93 astrometry were considered quite conservatively 
and include relative position measurement errors of +0713 
(1 pixel) as well as plate scale and detector orientation uncer- 
tainties of 070013 (1 % of plate scale) and +1.0°, respectively. 
While most of the pairs for which this comparison was possible 
shows no significant relative motion, a few cases retained our 
attention. The wide pair CVCha (~ 10" separation) is a clear 
outlier and the two components of this visual pair, CV Cha and 
CWCha, are most probably not physical. Hen 3-600 (TWA 3) 
is another binary which shows a significant relative motion and 
deserves further investigations (see AppendixlAl Another in- 
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Table 5. Probability of unrelated companions to the 
triple/quadruple systems as estimated from 2MASS. 





omp. 


[arcsec J 


P , j 

£ unrelated 


ESO Ha 283.... 


c 


1.36 10~ 3 


2.9 ur 2 




D 


1.70 10~ 3 


0.37 


PHty 30 


Q 


2.19 10~ 4 


8.8 ur 2 


VW Cha 


D 


1.27 10~ 5 


5.4 10~ 3 


CGHff 5/6 


B 


7.10 10- 6 


5.6 10-" 


C 


7.10 10- 6 


2.8 10~ 3 


N/ II 


C 


4.94 10- 6 


2.1 10~ 3 


Sz 30 


C 


1.45 10~ 5 


9.2 10~ 4 


UZ Tau 


B 


3.09 ur 7 


1.3 10 5 


.1 4872 


C 


5.25 10~ 6 


1.8 10~ 4 


UX Tau 


D 


7.10 10- 6 


1.6 10 4 


LkH« 336 


C 


1.67 10~ 5 


1.8 10- 3 


FV Tau 


c 


2.78 10~ 6 


1.3 10~ 3 


GG Tau 


c 


8.02 i<r 6 


2.6 10~ 3 


LkHa 262/263... 


D 


2.16 10 5 


1.6 10~ 3 


LkHa 346 


C 


1.54 10- 4 


5.9 10~ 4 




D 


3.70 10~ 5 


3.0 10- 3 


ROXs 43 


B 


i.88 ur 5 


5.3 10~ 6 




C 


2.16 10 s 


1.4 10" 4 


SR 24 


c 


1.36 10" 5 


1.1 10~ 3 


Sz 68 


B 


4.94 10~ 6 


2.4 10~ 7 




C 


1.27 ur 5 


3.1 10-" 



teresting point is that the seemingly non-hierachical systems 
LkHa346 (AB and AC pair ), Sz30 (AB and AC pair), and 
UX Tau (AB pair) does not show any significant relative mo- 
tion within the ~ 13 yrs time span. 

In some cases, it was possible to compare the inferred rel- 
ative motion with measured proper-motions (e.g. Ducourant 
et al. 12005b . In addition, among the four systems with re- 
solved proper-motions in the catalog of Ducourant et al. ( 2005 1, 
namely VWCha/Sz23, CVCha/CWCha, FV Tau/FV Tau, 
LkHa263/LkHa262, only CVCha/CWCha shows a signifi- 
cant (more than 3 sigma) relative proper-motion. Although the 
relative proper motion of the CV Cha system is mainly in dec- 
lination (also confirmed by Tycho-2 catalog), which is consis- 
tent with the relative motion we see here and would tend to 
favor the chance projection scenario, its amplitude would lead 
to a relative motion not exceeding ~U.'3 in 13 yrs. We there- 
fore suggest either an underestimate of the uncertainties on the 
proper-motions or an error in the astrometry reported in RZ93. 



4.5. Chance projections 

In order to infer what are the systems whose components are 
gravitationally bound and those which are only the result of 
chance projection, we used two approaches. 

The first one is a statistical approach which consists in esti- 
mating the probability that the companions we found are phys- 
ically bound to their primary. It implies to determine first the 
local surface density of background/foreground sources in each 
field. For that purpose we compiled the number of 2MASS ob- 
jects at least as bright in K-band as the candidate companion(s) 
in a 30x30 arcmin field surrounding each primary. This leads to 
the average surface density of objects brighter than this limiting 
magnitude, noted £ (K<K comp ). Assuming a random uniform 
distribution of unrelated objects across the field, the resulting 
probability P(X, 0) of at least an unrelated source to be located 




0.001 

0.1 1.0 10.0 

separation [arcsec] 



Fig. 4. Results of our survey in a plot of flux ratio or magnitude 
difference in H and K vs. pair separation. The circle symbol 
refers to pairs in triple systems while square symbols are for 
pairs in quadruple systems. The lines represent the complete- 
ness of our observations as derived from the sensitivity limit 
for undetected companions of the only binaries of the sample 
(see Sect.l4"2l. The lines are, from top to bottom, 99%, 90%, 
70%, 50%, 30%, 10% and 1% completeness. 

within a certain angular distance © from a particular target is 
given by : 

/>(£,©) = 1 -e^ 102 . (1) 

Since we did not attempt to distinguish between cloud mem- 
bers and background/foreground stars, the resulting probabil- 
ity is therefore also accounting for unrelated cloud members. 
The last column of Table[5]gives the resulting probabilities for 
a companion to be unrelated to the primary of a system. Most 
of the companions detected in our survey have probabilities to 
be projected unrelated stars well below the percent level. This 
means that they are very likely bound to their systems, although 
considering probabilities to individual sources is known to be 
prone to error (see e.g. Brandner et al. 2000 for a discussion). 
Three candidate companions (ESO Ha 283 C, ESO Ha 283 D, 
and PHa 30 C) show however a non-negligible probability of 
being the result of chance projections with probabilities of re- 
spectively 2.9%, 37% and 8.8%. 

The second approach is an attempt to determine the nature 
of the new or so far unconfirmed components through the use 
of both a color-color J-H/H-K diagram and a H-K/M^ color- 
magnitude diagram. While it was possible to composed the 
latter with all resolved components of the first data set, for 
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Fig. 5. Upper left : J-H/H-K color-color diagram for all sources with reliable 2MASS magnitudes, i.e. detected in all three bands, 
in the FOV of our targets. Blue small dot symbols denote sources detected in our survey while green triangle symbols are 
those non or marginally detected. Upper right : J-H/H-K color-color diagram for the new (PHa 30, CGHa 5/6, ESO Ha 283) and 
unconfirmed (LkHa336 and Sz41) systems. Filled squares represent component(s) included in the ~2 arcsec beam of 2MASS 
along with the main component, filled circles are possible companions linked to the latter by dotted lines. The loci of unreddened 
dwarf and giant stars are shown as the solid lines (Bessel & Brett 1988 ), as well as the locus of unreddened classical T Tauri stars 
(CTTS) as the dashed line (Meyer et al. U997t . A reddening vector for a standard reddening law is shown (E(J-H) = 0.11 Ay, 
E(H-K) = 0.065 Ay, Meyer et al 1997 1. Lower left : H-K/M^ color-magnitude diagram for all resolved components in the FOV 
of our targets for which we have a H-K color. As above, blue small dot symbols denote sources detected in our survey while 
green triangle symbols are those non-detected. The distribution of component is relative to the ZAMS, to the beginning of the 
giant branch (GB), and to a 1 Myr isochrone for PMS stars from 0.08 to 1 .4 M Q (with [M/H]=0, Y=0.275 and L mix = U p , Baraffe 
et al. U998> . Lower right: H-K/M^ color-magnitude diagram for all resolved components of the new and unconfirmed systems 
above. Absolute magnitudes M% are computed using the adopted distances shown in Tableland assuming 10% relative error. 
The reddening vector assumes A^ = 0.11 Ay. All data are in the CIT system. 



which we have H-K color measurements, the former is only 
composed of the well separated components for which 2MASS 
J, H and K magnitudes are available. Thus, 2MASS combined 
magnitudes and/or single-component magnitudes for the well- 
resolved (> 2 arcsec) companions are used for the color-color 
diagram. This means that through this latter we are probing the 
only well-separated candidate companions in our survey. This 
is nevertheless still useful since these candidate companions are 
usually those displaying the larger probabilities to be unrelated 
background sources (Table[5}. 



The five systems PHa 30, CGHa 5/6, ESO Ha 283, 
LkHa 336 and Sz 41 are plotted in the upper-right plot of Fig.|5] 
with the combined color of the brightest components as de- 
rived by 2MASS as filled squares. It turns out that PHa 30 C 
and Sz 41 C might well be giant background stars with very low 
line-of-sight extinction in view of this diagram. ESO Ha 283 D 
can also be interpreted as a background giant seen under 5- 
6 mag of visual extinction. However, since WTTS are also often 
located near the giant locus, as it can be noted in the upper-left 
plot of Fig.|5]where all systems have been plotted, one cannot 
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conclude on the sole basis of this diagram, but it provides fur- 
ther evidence for PHa30C and ESO Ha 283 D to be chance 
projections. 

The color-magnitude diagram presented in the lower-right 
plot of Fig. |5] is composed of all resolved components of the 
five sytems in our frames, while the color-magnitude diagram 
represented in the lower-left plot of the same figure stands for 
all systems with known H-K colour. Magnitudes for the sin- 
gle components are derived using H and K 2MASS combined 
magnitudes and relative photometry in [Fell] and Bry, respec- 
tively. This plot shows that, given the known distance to Sz 41, 
PHq-30 and ESO Ha 283 (respectively 160, 450 and 700 pc), 
the hypothesis that Sz41 C, PHa30C, ESO Ha 283 C and D 
are actually background giants implies distances of at least 
~480pc, ~4,5kpc and ~7kpc, respectively, which are consis- 
tent with giant absolute brightness. It has been actually shown 
that Sz 41 C is indeed a background giant (see AppendixlAT 

In summary, it is clear from these diagrams that some can- 
didate companions (especially the faintest objects) are likely 
extincted background stars. Others are consistent with T Tauri 
stars suffering differential extinction, or showing perhaps dif- 
ferent accretion disk properties and orientations. Although 
spectroscopy and common proper-motion are necessary in or- 
der to unambiguously identify any chance projection, we con- 
clude from the above analysis that PHa30C, ESO Ha 283 C 
and ESO Ha 283 D are consistent with being projected back- 
ground stars. As for Sz41 C, we will not consider further these 
companion candidates in our analysis. 

The compilation of all 2MASS sources with magnitudes 
in all JHK bands present in our FOV reveals some objects 
which are marginally detected (55 5-sigma) or undetected in our 
frames. These sources are reported in the plots of Fig.[5]with 
triangle symbols. These 7 sources are in the FOV of CV Cha 
(2 objects), Sz62, PHa 14, PHa30, CGHa5/6 and GGTau. 
These sources are mostly consistent with background main- 
sequence stars or giants in view of their locations on these di- 
agrams. In addition, it is most probable that they would have 
been identified as cloud members by previous H a surveys, even 
if we cannot totally exclude that some could still be very-low 
mass (550. 3Mq) young (SSlOMyr) non-accreting stars, there- 
fore difficult to be identified by such surveys. A more complete 
census of the low-mass content of clouds as the one provided 
by e.g. Luhman (2004i for Cha I, in association with resolved 
X-ray studies would certainly shed light on the nature of these 
sources. Finally, it is noteworthy to point out that some systems 
present unsual bluer colors and are located to the right of the 
main-sequence in both diagrams. These sources corresponds to 
J 4872 ABCD (only resolved in the color-magnitude diagram), 
UXTauBC (idem) and Sz60A. While further investigations 
would be necessary, especially resolved spectroscopy and JHK 
photometry, we suspect that the proximity with the other com- 
ponents in those systems (with separations S5 3'.'5) might have 
biased the 2MASS magnitudes. 




Log(Te,f) 

Fig. 6. Locus of new and poorly studied companions in the 
(T e ff , Mh) plane together with the PMS evolutionary tracks and 
isochrones of Baraffe et al. ( 1998 ). The components are placed 
in the diagram using their Mh and assuming coevality with the 
closest component of the system. 

4.6. Nature of the companions 

We derived an estimate of the spectral types and masses of the 
new, presumably physical, companions as well as those that are 
poorly studied (and lack of spectral type and/or mass estimates) 
from their H-band magnitude and the assumption of coevality 
with the other stars in the systems. These values are reported in 
Tabled 

With the H-band magnitude from Tableland spectral type 
from the literature we placed the known and studied mem- 
bers of each system on the (T e ff, Mh) plane of the H-R di- 
agram with PMS evolutionary tracks and isochrones overlaid 
and derived an age for the systems. The values of T e ff for 
these known and studied components of the systems are de- 
rived from the system's spectral type quoted in literature (re- 
ported in Table[Q and the temperature scale from Sherry et 
al. l|2)04). When T e ff values are quoted in literature they are 
within half a subclass the one derived with the adopted temper- 
ature scale. Uncertainties of typically one subclass in spectral 
type are assumed, apart for some cases for which uncertainties 
could be much larger (e.g. CGH„ 5/6 B). Uncertainties in abso- 
lute H-band magnitude assume a distance uncertainty of 10%, 
and values of Ay from the literature were used to deredden 
the H-band magnitudes. In Fig.[6]we show the only new multi- 
ple system components as well as poorly studied components 
plotted on the H-R diagram together with the PMS tracks and 
isochrones of Baraffe et al. ( 1998) with a=1.9 for M>0.6M Q . 
A similar diagram not shown here was formed for Sz 30 C for 
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which only K-magnitudes are available. Uncertainties in the 
inferred masses and spectral types of these new and/or poorly 
studied components are set by changing the absolute H-band 
magnitude by its uncertainty and searching for the intersection 
with the isochrone corresponding to the inferred age. An ex- 
ample of the procedure in the case of LkHa 346 is presented in 
AppendixlAl Due to the fact that H-band flux might not be com- 
pletely photospheric (Greene & Meyer 1995 1, the spectral types 
and a fortiori the masses derived here should only be consid- 
ered as crude estimates. This is especially the case for Sz30C 
for which we had to use K-band fluxes. We would therefore 
expect an underestimate of the ages of the primaries which re- 
flects on earlier spectral-types and overestimated masses for the 
assumed coeval companions. Future follow-up spectroscopy of 
these objects will give much more accurate T e ff values allow- 
ing to test coevality of these system components. In this re- 
spect, if the two new components of the LkHa 346 system are 
confirmed to be physical, then this system would be an im- 
portant test case for the PMS evolutionary models, spanning 
a relatively large range of masses and possibly extending to 
the substellar regime, in the same fashion as GG Tau (White et 
al.fT999l. 



4.7. Multiplicity statistics 

Among the 58 wide binaries surveyed, five are Herbig 
Ae/Be binary stars (HD 76534, Herschel 4636, PHa 30, Sz 120, 
Elias 2-49) and one is likely to be a foreground (older) object 
(Sz 15, see AppendixlAl. Excluding these systems from our 
statistics, we have 52 T Tauri star binary systems with mostly K 
and M-type primary, apart COOri (F8e), Sz 19 (G2), CVCha 
(G8), Elias 2-30 (G2.5) and ROXs 43 (GO). Taking into account 
the edge-on disk LkHa 262/263 C we end up with 38 binaries, 
7 triples (CGHa 5/6, LkHa 336, ROXs 43, SR 24, Sz 30, Sz 68, 
UZ Tau) and 7 quadruples (LkHa 262/263, LkHa 346, FV Tau, 
GG Tau, J 4872, UX Tau, VW Cha). We did not attempt to cor- 
rect for incompleteness (see discussion in Sect. l4.2t . In order 
to characterize the multiplicity, we here define two quantities. 
First, a quantity we call degree of multiplicity per wide binary 
(or a multiplicity frequency per wide binary, MF/wB) : 



MF/wB = 



T + Q + 



wB + T + Q + 



(2) 



where wB represents the number of wide binaries (with pro- 
jected component separations typically 075, see Sect.[2), T 
the number of triples and Q the number of quadruples. This 
quantity here equals 26.9+7.2%. The errors have been esti- 
mated by taking the square root of the number of triples and 
quadruples. 

Second, a so-called companion star frequency per wide bi- 
nary (CSF/wB) which is defined as : 



CSF/wB = 



27 + 3(2 + .. 
wB + T + Q + 



(3) 



Our study gives 0.67 +0. 1 1 for this quantity. Note that the lat- 
ter quantity is sensitive to a decrease of high-order multiplicity 
by, e.g., dynamical decay while the former quantity is not. In 



general, while MF/wB is only a measure of how common high- 
order systems are, CSF/wB also give us information about the 
relative fractions of high-order multiplicity via the number of 
companions. 

In the rest of the study, we will consider a distance-limited 
sample in order to ensure a similar range of linear projected 
separations probed, thus providing meaningful statistical re- 
sults. We limit ourselves to only the wide T Tauri binaries of 
the sample at distance 140-190pc (41 systems). This includes 
only the systems for which the separation between compan- 
ions are in the range 10/14 AU - 1700/2300 AU, i.e. with 
projected separations between 0707 and 12". In that case, 
MF/wB=26.8±8.1% and CSF/wB=0.68 ±0.13 (30 binaries, 5 
triples, 6 quadruples). 



5. Discussion 

5.1. Multiplicity as a function of mass 

A breakdown of the distance-limited sample by primary mass, 
i.e. by spectral type of the primary, leads to MF/wB=38 +14% 
and CSF/wB =1.1 +0.2 for K types (13 binaries, 2 triples, 6 
quadruples), MF/wB=13+9% and CSF/wB=0.3 ±0.1 for M 
types (13 binaries, 2 triples). The trend is consistent with the 
known decrease of multiplicity with decreasing primary mass 
(Sterzik & Durisen l2003l . 

5.2. Comparison between clouds 

The question that arises naturally is the one about the multiplic- 
ity frequency in different clouds, although the latter is of lower 
statistical significance than that of the total sample. The cloud 
with the highest value is Tau-Aur (5 triples-quadruples/1 1 
wide binaries, MF/wB =45 ±20%, CSF/wB=1.3 ±0.3), fol- 
lowed by Chal (2/10, MF/wB=20±14%, CSF/wB=0.5 ±0.2) 
and Ophiuchus (2/10, MF/wB=20±14%, CSF/wB=0.4 ±0.2). 
Some star formation models suggest that cloud intrinsic prop- 
erties might imprint different binary or multiplicity properties 
(Sterzik & Durisen [T9"9"8l Sterzik et al. l2003> . but no model is 
able to quantify the differences observed. 

5.3. Comparison with previous multiplicity surveys 

We compare our result with the proportion of 
triples/quadruples found in previous NIR multiplicity surveys 
which have focussed on the brightest sources (generally 
K $ 10) with broadly similar separation range, sensitivity 4 , and 
statistical uncertainties. These are the studies by Leinert et al. 
(Q3533J! and Kohler & Leinert ( fTWST l in Tau-Aur, Ghez et al. 
H1997a> in Chamaeleon, Lupus and CrA, Kohler et al. J2000b> 
in the Scorpius-Centaurus OB association, Kohler (2001 1 
in Chamaeleon, and Ratzka et al. (2005 1 in Ophiuchus. We 
primarily base this comparison on the multiplicity frequency 
per wide binary (MF/wB), as defined above (Eq|2|l, and 
the companion star frequency per wide binary (CSF/wB, 



4 However, our AO survey is probably slightly more sensitive than 
previous speckle surveys. 
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Table 6. Comparison of the multiplicity frequency per binary (MF/wB and MF/B) and the companion star frequency per binary 
(CSF/wB and CSF/B) of our work with those derived from previous multiplicity surveys among T Tauri stars in a similar 
separation range ~14-1700 ALT", and with recent numerical simulations (bottom part). Also indicated are the number of wide 
(separation > Of! 5) binaries (wB), binaries (B), triples (T), quadruples (Q) and higher multiples (>4). 



Reference 


cloud 


wB 


B 


T 


Q >4 


MF/wB 


CSF/wB 


MF/B 


CSF/B 


This work 


several 


30 


42.3 


5 


6 


26.8 ±8.1% 


0.68±0.13 


20.6 ±6.2% 


0.53±0.10 










Other Multiplicity Survey: 










L93 


Tau-Aur 


24 


39 


3 


2 


17.2 ±7.7% 


0.41 ±0.12 


11.4*5.1% 


0.27 ±0.08 


KL98 


Tau-Aur 


21 


30 


6 


1 


25.0 ±9.4% 


0.54 ±0.14 


18.9 ±7.2% 


0.41 ±0.11 


G97 


Cha/Lup/CrA 


19 


21 


2 




9.5 ±6.7% 


0.19±0.10 


8.7 ±6.1% 


0.17 ±0.09 


K00 


Sco-Cen 


24 


37 


6 




20.0 ±8.2% 


0.40 ±0.12 


14.0 ±5.7% 


0.28 ±0.08 


KOI 


Cha 


5 


16 


1 




16.7 ±16.7% 


0.33 ±0.24 


5.9 ±5.9% 


0.12±0.08 


R05 


Oph 


22 


45 


5 




18.5 ±8.3% 


0.37±0.12 


10.0 ±4.5% 


0.20 ± 0.06 


Tot(Other surveys) 




115 


188 


23 


3 


18.4 ±3.6% 


0.39 ±0.05 


12.2 ±2.4% 


0.26 ±0.03 


Tot(All surveys) 




132 


217.3 


23 


7 


18.5 ±3.4% 


0.41 ±0.05 


12.1 ±2.2% 


0.27 ±0.03 










Numerical Simulations 










SD03 






0.21 


0.08 


0.03 






34% 


0.78 


D03 






0.154 


0.108 


0.045 






49.8% 


1.14 


D04 






11 


2 


2 3 






38.9 ±14.7% 


1.28 ±0.27 


G04 






7 


6 


2 1 






56.3 ±18.8% 


1.38 ±0.29 



a : This range is assuming that all systems are at 140 pc. Additionally the upper limit is actually set by confusion with background/foreground stars, and varies from survey to 
survey. The surveys of Leinert et al. 1993 i, Kohler & Leine rt 119981 . and Ghez et al. I1997al have set their upper limit to 12" (~ 1700 AU at 140 pc), while the surveys of 
Kohler et al. I2000bl . Kohler |2001land Ratzka et al. 120051 h ave up per limits of 6". 

References : L93 : Leinert et a UlTOl KL98 : Kohler & Lein ert fiWl G97 : Ghez et al. H997al K00 : Kohler et al. 2000b K01 : Kohle nTOTl R05 : Ratzka et al.[5005; SD03 : 
Sterzik & Durisen l2003l D03 : Delgado-Donate et al. 120031 D04 : Delgado-Donate et al. l2004l G04 : Goodwin et al. 120041 



Eq|3J- Since those surveys used samples drawn from lists of 
apparently single stars, and not wide binaries like our, they 
include all resolved binaries with separations down to the 
resolution limit. They therefore give multiplicity fractions 
per binary, i.e. MF/B and CSF/B, where B includes binaries 
with separations S5075. We thus have to count in these surveys 
only binaries with separations larger than 075 (i.e. ~70 AU at 
140 pc). This will give us an estimate of MF/wB and CSF/wB 
for these surveys, and allow a direct comparison with our 
own survey results. Note that in all these surveys we consider 
a companion number uncorrected for either completeness 
or chance projection, as these corrections are difficult to 
assess properly and moreover roughly compensating each 
other. While completeness depends on the sensitivity of these 
surveys which is rather uniform, chance projection is most 
of the time treated by assuming a 'safe' maximum separation 
which depends on the background/foreground surface density 
of stars. On our side, we would have to slightly restrict the 
separation range to the resolution achieved by those surveys 
(i.e. typically -071-12" which is -14-1700AU at 140pc), 
but this is not significant since only one companion was 
detected below this limit (SR24C at separation 0708). Tabled 
summarizes the result of this comparison. It shows that in first 
order our newly derived multiplicity agrees with the previous 
surveys, within the uncertainties. 

In a second step, we will attempt to estimate the fraction of 
binaries which were not included in our survey, i.e. those with 
a separation range ~07 1-075. This will allow us to convert our 
MF/wB to MF/B (as well as CSF/wB to CSF/B) and to perform 
a direct comparison with the results from the other surveys. For 
that purpose, we derived the fraction of binaries of those sur- 
veys with separations between ~071 and 075. We only consider 
those surveys for which the maximum separation was set to 
12", and find a fraction of close binaries a=28.9 ±5.7%. The 



difference with the fraction derived from all surveys, including 
those with upper separations of 6", is not statistically differ- 
ent since in that case Q'=38.8 ±4.5%. The corrected number of 
binaries with separations >071 for our survey is then B=wB/(l- 
a)=30/(l-0.29)=42.3. Thus, MF/B = 20.6 +6.2% and CSF/B = 
0.53 ±0.10 leading to the same conclusion as above. 

Noticeably, the CSF/B (and CSF/wB) is significantly 
higher in our survey than for all the other surveys. If we av- 
erage the result for all the other surveys the difference is at the 
2<x level. This overabundance of companions is a consequence 
of a higher fraction of quadruples over triples in our survey 
with respect to all previous surveys. The fraction of systems 
with multiplicity higher or equal than four to the number of 
systems with multiplicity higher or equal than three (called f\ 
in the Sect. l5.5t is 55.0 ±22.2% for our survey, while it is only 
13.0 ±7.5% on average for all other surveys. 

5.4. Comparison with theory 

There is a probable overabundance of high-order multiples pro- 
duced by the current simulations of star formation with respect 
to current observations. Direct comparison is not possible since 
theoretical multiplicity frequencies include binaries with sep- 
arations down to typically ~ 3-5 AU. Wider high-order com- 
panions with separations ^ 2000 AU are quite rare, so that the 
upper bound of the separation range should be less of a con- 
cern. We will assume here that the corrections to be applied 
in order to obtain MF/B in the same separation range as the 
one probed here are minor, i.e. we will neglect the correction 
for systems in the separation range ~ 3-5 AU-10AU. In some 
cases, this might not be the case. For example in the simula- 
tions of Goodwin et al. (20041 the fraction of close binaries 
with separations ^ 10 AU would be up to half of all binaries 
formed (S. Goodwin 2006, private communication). We there- 
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fore have to keep in mind that the theoretically derived MF/B 
should be considered as lower-limits. In the following, we sum- 
marize the theoretical studies used for the comparison. 

Sterzik & Durisen ( 2003 1 performed few -body cluster de- 
cay simulations. Although that study neglected the effect of 
remnant molecular gas and disk accretion and treated only the 
process of dynamical evolution of young small N-body clus- 
ters, it yields highly significant and robust statistics since a 
large number of realizations (10 000) has been computed. A de- 
gree of multiplicity of MF/B=34% (CSF/B = 0.78) was found. 
Delgado-Donate et al. (2003 1 modeled the dynamical decay 
of a large number (a hundred) of small-N (N=5) star-forming 
clusters including the effects of competitive accretion and dy- 
namical evolution through Smoothed Particles Hydrodynamics 
(SPH) simulations with a ~1AU spatial resolution, and found 
a rather high multiplicity frequency close to MF/B =50%. A 
similar high frequency of multiple systems was the outcome 
of two other recent and more sophisticated SPH simulations. 
Delgado-Donate et al. (2004 1 simulated the fragmentation of 
10 small-scale turbulent molecular clouds and their subsequent 
dynamical evolution, including this time the effect of accretion 
disks into the evolution of multiples. Goodwin et al. (2004i fol- 
lowed the collapse and fragmentation of 20 dense star-forming 
cores with a low-level of turbulence. In both cases a high fre- 
quency of high-order multiples was obtained (Table[6j. 

5.5. Comparison with Main-Sequence multiplicity 

Most of our targets will evolve towards G and K-types on 
the Main-Sequence, therefore we will compare our results 
with the field star surveys focussing on those spectral types. 
Duquennoy & Mayor (1991 hereafter DM91) gives an esti- 
mate of the high-order multiplicity among solar-type Main- 
Sequence dwarfs. However selection effects may apply, as al- 
ready realized in DM91, and more recent surveys and compila- 
tions (e.g. Tokovinin & Smekhov 2002 ) argue for a larger frac- 
tion of high-order multiples among MS dwarfs. In the follow- 
ing, we will compare the multiplicity fractions from the MSC 
catalog (Tokovinin 119971 with ours. MSC gives their results 
using a multiplicity fraction f„ defined as (Batten fT973> : 



where N n is the number of systems of multiplicity at least n. 
Batten (I1973i estimated f n to be of the order of 25%, inde- 
pendently of the order of multiplicity. Tokovinin (2006| es- 
timated from MSC a fraction of triple systems (or higher) 
relative to binaries (or higher) f 3=19+4% and a fraction of 
(at least) quadruples to (at least) triples /4=22±2%. A simi- 
lar value of /3=20+6% was found by Tokovinin & Smekhov 
( 2002) in a sample of dwarfs with spectral types later than F5 
extracted from MSC. This latter value is very similar to the for- 
mer one, suggesting that the high-mass evolved stars included 
in MSC should have little effects on the multiplicity fractions. 
In comparison, a value of /3 = 11 +4% can be estimated from 
the sample of DM91. In our case, it happens that /3=MF/B, 
so that / 3 =20.6±6.2% and / 4 =55.0 ±22.2%. While the frac- 
tion of triples to binaries is in agreement with what is found for 



the main-sequence, there is a marginally significant overabun- 
dance of quadruples with respect to triples in our survey when 
compared to MSC values. A clear caveat of this comparison is 
that multiplicity fractions for main-sequence stars are derived 
for much larger linear separation range than for our survey. We 
attempt to address this issue in the next section. 

5.6. Total multiplicity 

We searched in the literature for spectroscopic binaries (SBs) 
included in our sample. Systematic searches for spectrsocopic 
companions of T Tauri stars comprise the radial-velocity sur- 
veys of Mathieu et al. ( 119891 . Reipurth et al. (2002), Torres 
et al. ( I2003I I. and Melo 12003t . According to them, there are 8 
SBs, of which 5 need confirmation, in our sample (see TableQ. 
Since not all the components of our sample have been searched 
for spectroscopic companions, this number is likely to be a 
lower limit. We can therefore tentatively estimate the increase 
of multiplicity that would result if we could include SBs by 
taking into account all those potential 8 SBs. This would lead 
to multiplicities MF/wB=39.0 +9.8% and CSF/wB=1.02 +0.16 
for companions with separation down to a few tens of AU 
(25 binaries, 8 triples, 6 quadruples, 2 quintuples) 5 . This 
also corresponds to multiplicity fractions (Eq.@} /3=39 ±10%, 
which is an upper limit here (see below), /4=50±18%, and 
/s=25 ±18%, i.e. the ratio of quadruples to triples is still higher 
than that of the main-sequence (f$ is within the uncertainties of 
that of MSC which is 20 ±4%). 

We try to better estimate by tentatively correcting 
for the binaries with separations lower than 075 as before, 
but not for those systems with separation between a few 
tenths of AU to ~10AU, i.e. almost pure SBs, as these 
systems might not be so ubiquitous (Sterzik et al.|2005 
Tokovinin et al. l2006l . This gives MF/B=/ 3 =31.3 ±7.8% and 
CSF/B=0.82±0.13. Therfore, the ratio of triples to binaries 
appears also higher for young systems than for the main- 
sequence, although at a marginally significant statistical level 
and with the caution that this result might be in part biased 
by the assumption we made about the paucity of tight binaries. 
Noticeably, a relatively good match is found with the numerical 
models of dynamical evolution of young small N-body clusters 
(see Sect.Eland TablegJ. 

5.7. Hierarchical configurations 

Among the quadruple systems, three systems are composed of 
two pairs separated by a distance larger in projection than the 
pair separations, and four systems are composed of a pair with 
a wide companion and an even wider companion. An almost 
equal number of systems between these two hierarchical con- 
figurations is also what is found among Main-Sequence sys- 
tems in the MSC catalog (Tokovinin 20(H). 



5 This includes the 07016 companion to ROXs43C (Simon et 

ai. rn^5i 



16 



S. Correia et al.: A VLT/NACO Survey for Triple and Quadruple PMS Systems 



Table 7. Systems dynamical stability. Application of the sta- 
bility criteria from Eggleton & Kiseleva (1995) (their Y™", 

cf Eq. 2) and from Mardling & Aarseth (2001) (R™ = ^~ 
their Eq. 90) for a given partition (l+2)+3, q,„=mi/m2 and 
(\om-(fn\ + m2)/ni3. R mes is the ratio of the outer to the inner 
projected separation of the triple partition considered". 



System Partition 6 q,„ q„,/ Y™""' R™.™<" R, 



VWCha (B+Q+A 1.14 0.75 6.16 3.92 5.8 

CGHo-5/6 (A+B)+C 1.84 0.93 5.89 3.75 23.4 

Sz30 (A+B)+C 1.60 1.34 5.45 3.49 3.1 

UZTau (B+Q+A 1.08 0.60 5.25 3.38 9.8 

J 4872 (A+B)+CD 1.10 1.39 5.43 3.48 19.3 

UXTau (A+D)+BC 6.81 2.40 4.81 3.22 2.2 

LkHo-336 (A+C)+B 2.83 1.33 5.46 3.5 1.41 

FVTau (C+D)+AB 3.15 0.27 7.87 5.20 17.4 

GGTau (C+D)+AB 2.85 0.11 9.80 .J 6.92 

LkHff 262/263... (A+B)+C 0.95 1.17 5.62 3.58 9.9 

LkHa-346 (A+B)+CD 3.11 0.42 7.08 4.57 5.3 

ROXs43 (A+B)+C 1.13 1.16 5.63 3.59 13.4 

SR24 (B+Q+A 1.79 0.63 6.40 4.09 62.6 

Sz68 (A+B)+C 3.77 25.3 3.26 2.84 22.3 



a : In case of ambiguity the lowest value is adopted. 

b : When several configurations exist, the more unstable is reported. 

c : When several mass estimates exist, these columns report the values corresponding to the 

greater values of Y^ and R™ f ". 

d : Not applicable since m^l(m\ + ni2) > 5. 

5.8. System dynamical stability 

In order to quantitatively investigate the system dynamical sta- 
bility, we applied to our sample the semi-empirical criteria 
developed by Eggleton & Kiseleva (1995) and Mardling & 
Aarseth (2001). The result is shown in Table0 where R mes 
should be at least larger than Y™" or R™" to ensure stability. 
Assumptions are that the periastron of the outer orbit as well 
as the semi-major axis and the apastron of the inner orbit are 
approximated by their projected separations, i.e. the excentric- 
ity of the outer orbit is assumed to be zero and the system is 
supposed to be seen face-on. The component masses used are 
taken from the literature when available and quoted in Tabled 
or roughly estimated by us using flux ratios and the assump- 
tion of coevality (Sect. l4.6t . Under the above assumptions, 3 
systems (LkHa336, Sz30 and UXTau) appear to be unstable 
by both criteria while three others (LkHa 346, VWCha and 
GGTau) appear unstable by the Eggleton & Kiseleva (1995) 
criterion only. In that respect, it is interesting to point out that 
recently, Beust & Dutrey ( 2006 ) have shown through a detailed 
analysis of the dynamics of the GG Tau system that only some 
special orbital configuration would allow GG Tau CD to be sta- 
ble, although this result depends strongly on the adopted size 
of the circumbinary disk around GG Tau AB. 

Even without the values reported in Table0 one can im- 
mediately note from Fig. |2 that the three systems LkHci'336, 
Sz30 and UXTau are apparently non-hierarchical. This cor- 
responds to a frequency 3/14=21 + 12%. Note that Kohler & 
Leinert ( 1998) found 1 apparent non-hierarchical triple among 
7 triples, while Kohler et al J2000bt and Ratzka et al. (2005) 
foundl/6 and 1/5, respectively. This yields a fraction of appar- 
ent non-hierarchical triples among triples of ~15-20%, and in 
combination with our survey, a total fraction of 6/32=19±8%. 

Ambartsumian (1954) provided a detailed statistical anal- 
ysis of how many Trapezium systems are in fact pseudo- 



Trapezium systems, i.e. projected hierarchical systems. With 
the assumption of circular orbits, the probability for a hierar- 
chical triple system to be seen in a pseudo-Trapezium configu- 
ration is (his Eq. 10) : 

u = n 1 1 ■ (5) 

Igki - lgk 

where ko < AC/AB < k\, AC and AB are the maximum and 
minimum projected separations, respectively, ko and k\ their ra- 
tios in the case of a Trapezium and hierarchical configuration, 
respectively. In our multiple system sample, we have ko ~ 2 and 
k\ ~ 10, thus U=0.05. Following Ambartsumian ( 1954 1, we in- 
clude the quadruple systems by considering an equal number of 
the two hierarchical configurations discussed in Sect. l5.7l and a 
fraction of quadruples among systems with multiplicity higher 
than 2 of one half. This leads to a probability P=0.07 that appar- 
ent non-hierarchical systems are in fact projected hierarchical 
systems and this is in statistical agreement with the value we 
found. 

However, if we suppose that the relative orbital inclina- 
tions could be rather small, as the study of Sterzik & Tokovinin 

(2002) among field triple systems seems to indicate, then the 
fraction of projected hierarchical systems among apparent non- 
hierarchical systems would be much smaller and 1 or 2 of our 3 
systems could be really non-hierarchical. These systems would 
then be either very young, or in re-arrangement after a recent 
close triple approach following or not the ejection of a compan- 
ion. In all cases, these systems are not expected to be long-term 
stable and may later decay. For example, Delgado-Donate et al. 

(2003 ) found that the non-hierarchical systems formed in their 
simulations decay in less than ~0.5Myr. UXTau is the only 
systems with an estimated age for A, D and the combined BC 
system (see Tabled, and it seems not consistent with a very 
young ($1 Myr) system. On the other hand, this could indeed 
be the case for LkHa 336. 

5.9. Disk evolution 

The study of disk evolution in young multiple system allows to 
investigate samples of stars of various masses that are a priori 
coeval. It is however necessary to first understand the influence 
of multiplicity on disk evolution, i.e. in which conditions the 
disks can be considered to evolve in isolation in multiple sys- 
tems. This has been the subject of several early studies among 
T Tauri binaries (see Monin et al. 2006b and references therein 
for a recent review). Only a few attempts to investigate disk 
evolution of higher order multiple systems have been achieved 
recently (White et al. 2002 1. This is however of particular inter- 
est since both close and wider pairs are usually present in such 
hierarchical systems and all components of a system are sup- 
posed to be coeval. In the following, we analyze the frequency 
of disks in our sample of triples/quadruples through the use of 
various disk and/or accretion diagnostics and compare to what 
is known for isolated binaries. 

We performed a compilation of T Tauri types for both indi- 
vidual components and close pairs of the triple/quadruple sys- 
tems from the available literature (Table|SJ. Two criteria were 
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Table 8. Known T Tauri types (CTTS or WTTS) for all single components or otherwise pairs of the triple/quadruple systems. 
Also included are spectral types, Ay, masses, and ages of the components from the literature, as well as our own estimates for 
the new and poorly studied components when possible (see Sect. l4.6l l. 
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used in order to assign these types : first, the equivalent width of 
the Ha line with a threshold between WTTS and CTTS at about 
10 A, depending on the spectral type (Martin 1998 Hartigan 
& Kenyon2003i; second, the A(K-L) color excess which was 



shown to correlate with both strong [OI] and Ha emission lines 
(Edwards et al. U993t . In the latter case A(K-L)>0.4 indicates 
the presence of an optically thick accretion disk (Edwards et 
al. 1993 1. The A(K-L) color excess was computed from K-L 
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measurements and Ay available in literature (for which the ref- 
erence is reported in Table[3}, spectral types either from the 
literature or estimated by us and intrinsic colors of Kenyon & 
Hartmann dl995> . For some systems individual component T 
Tauri types are missing and only the global type of the usually 
tightest pair is available. However, one can consider that a tight 
pair presenting spectral WTTS properties is in fact composed 
of two WTTS, while no conclusion is possible in the case of a 
tight pair for which the composite spectrum shows CTTS prop- 
erties. We recently obtained spatially-resolved spectroscopy for 
some of these tight systems and the results will be presented 
elsewhere. 

While the number of systems in Tabled does not allow us 
to draw any firm conclusions, some trends can nevertheless be 
identified. Among the 9 systems with CTTS/WTTS informa- 
tion known and/or deduced for each component, one is presum- 
ably composed of only WTTS (J 4872). Since this system does 
not contain information about differential disk evolution, we 
will focus on the remaining 8 systems (CGHa5/6, VWCha, 
UXTau, UZTau, FVTau, GGTau, LkHa346, SR24). Half of 
these 6 are systems of mixed type (i.e. at least one component 
with a different type). The fraction of mixed systems is even 
higher if one includes the 4 systems that are already mixed even 
without the full knowledge of all components' T Tauri types 
(i.e. Sz 30, LkHa 336, ROXs 43, Sz 68). We end up with a frac- 
tion of mixed systems of ~ 65 % (8/12), similar or even higher 
to what is found among binaries (~ 40 %, Monin et al. 2006b I. 
The much higher fraction of Taurus systems among non-mixed 
systems (3/4, compared to 1/8 for mixed systems) supports the 
suggestion of a regional dependance of the fraction of mixed 
systems, or at least a lower fraction in Taurus as in the case of 
binaries (Monin et al. 2006b ). On the other hand, we do not find 
any difference in the distribution of estimated age in mixed- 
and non-mixed systems which could explain the existence of 
these two populations, i.e. the age of non-mixed systems is not 
significantly younger than that of mixed systems. 

Disk diagnostics in all close pairs (separation $ 100AU) 
are consistent with what is expected from tidal truncation the- 
ory (Armitage et al. fT^t . All close CTTS pairs (GGTau AB, 
UZTau BC, SR24BC, FVTau AB) exhibit estimated mass ra- 
tios qS0.5. SR24BC with q~0.5 has the lowest value, and 
this particular case may need disk replenishment from a com- 
mon envelope with material falling preferentially onto the sec- 
ondary's disk (Bate & Bonnell 1997 1 in order to prevent the 
more truncated of the secondary disk to be dissipated prior 
that of the primary. In general, all four close CTTS pairs may 
well need disk replenishment in order for the highly trun- 
cated disks to sustain active accretion. While both GG Tau AB 
and SR24BC are known to be surrounded by circumbi- 
nary disks for replenishment (Dutrey et al. l 19941 Andrews & 
Williams 2005), the cases of FVTau AB and UZTauBC seem 
less obvious since no circumbinary disk has been detected so 
far. A massive circumbinary disk is known to be surrounding 
the spectroscopic pair UZTau A located some 500 AU from 
UZ Tau BC but it is very unlikely responsible for the resplen- 
ishment of the BC pair. It is interesting to note that only gas 

6 FV Tau, as discussed below in the text, is not a mixed system. 



emission in the CO line was detected around the pair SR 24 BC 
(Andrews & Williams 2005 1 which suggests that other cir- 
cumbinary disks or common reservoir could have been missed 
by previous continuum millimeter surveys. The case of the 
FVTau CD pair deserves special comment since it seems to 
be, together with VW Cha A-BC, a close pair of mixed type. 
While the latter is consistent with exhausted highly truncated 
secondary disks, FVTau CD shows a WTTS-CTTS type de- 
spite an estimated q~ 0.25 which would imply a more tidally 
truncated secondary disk. However, one should note that this 
apparent contradiction can be solved if one considers that com- 
ponent C is in fact a passive disk (i.e. a non-accreting disk, 
e.g. McCabe et al.2006) and component D an IRC candidate 
(Woitas et al. 120011 see also McCabe et al. l2006t . Therefore, 
while the origin of the IRC phenomenon is still under debate, 
we can speculate that the latter is likely to be more massive, i.e. 
the primary. This implies that the close pair FV Tau CD which 
has two disks would have a star mass ratio presumably closer 
to unity, hence in agreement with synchronized disk evolution. 

The fact that 4 mixed systems include both close WTTS 
pairs and well-separated components with disks seems to sug- 
gest that the highly truncated disks of these close pairs were 
not fed by a circumbinary envelope and hence became rapidily 
exhausted. Alternatively a reservoir could have originally ex- 
isted in such multiple systems and be subsequently disrupted 
by dynamical interactions with the other component(s) of these 
systems. 

6. Summary 

In this paper, we reported on our survey for high-order multi- 
plicity among wide visual Pre-Main Sequence (PMS) binaries 
conducted with NACO at the VLT. The main conclusions of 
our study are summarized as follows : 

(1) Among the 58 PMS wide binaries surveyed, which com- 
prises 52 T Tauri systems from various star-forming regions, 
we found 7 triple systems (2 new) and 7 quadruple systems 
(1 new). The new close companions are most likely physically 
bound based on their probability of chance projection and on 
their position on a color-color diagram. Some systems might 
still be not physical and future spectroscopy and common- 
proper motion measurement will be able to give an answer 
soon. The corresponding degree of multiplicity among wide bi- 
naries is MF/wB=26.9 +7.2% in the projected separation range 
0707-12", with the largest contribution from the Taurus cloud. 
Considering a restricted sample composed of systems at dis- 
tance 140-190pc, we obtained only a slightly different value 
of MF/wB=26.8+8.1%, in the separation range 10/14 AU - 
1700/2300 AU (30 binaries, 5 triples, 6 quadruples). 

(2) A general trend of decreasing multiplicity with lower pri- 
mary mass is found, consistent with the known decrease of 
multiplicity with primary mass mainly found among binaries 
(Sterzik & Duri sen 120031 . 

(3) Comparison with previous multiplicity surveys focusing 
on the brightest sources with similar resolution and (perhaps 
slightly lower) sensitivity shows that there is a good agree- 
ment between our newly derived multiplicity frequency per 
wide binary and those derived from these surveys, in the com- 
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mon separation range ~14-1700 AU, although a significant 
overabundance of quadruple systems compared to triple sys- 
tems is apparent. Considering that 4 out of the 6 quadru- 
ples of our distance-limited sample are from Taurus and that 
sources of that region have a frequency of quadruples to triples 
/4=39+17% when combining the results from previous sur- 
veys with ours, indicates that this excess of quadruples is 
mainly due to Taurus. 

(4) Tentatively including the spectroscopic pairs to our re- 
stricted sample and comparing the multiplicity fractions to 
those measured from solar-type main-sequence stars in the so- 
lar neighborhood leads to the conclusion that both the ratio of 
triples to binaries (/j) and the ratio of quadruples to triples (J4) 
seems to be in excess among young stars. While the former is 
not statistically significant, and depends on the assumption of 
paucity of spectroscopic binaries with no third component, the 
latter may be statistically significant. 

(5) Our multiplicity frequency per binary seems lower than cur- 
rent predictions from numerical simulations of multiple star 
formation, especially SPH simulations. Noticeably, however, 
a relatively good match is found with the numerical models of 
the dynamical evolution of young small N-body clusters. 

(6) Some systems might not be dynamically stable according 
to the criteria of Eggleton & Kiseleva ( 1995 1 and Mardling & 
Aarseth (2001 1. However, the three apparent non-hierarchical 
triple systems could be real if one assumes some correlation in 
the orbital inclinations. In one case, the age of the system could 
be young enough to explain this configuration. 

(7) We performed a compilation of T Tauri types for both in- 
dividual components and pairs of the triple/quadruple systems 
from the available literature in order to study the relative disk 
evolutions in these systems. With the caution of small num- 
ber statistics we identified the following trends. The fraction of 
mixed systems is high, up to ~ 65 %, similar to what is found 
for binaries. Likewise, the fraction of mixed systems in Taurus 
is low. Disk replenishment by a circumbinary envelope seems 
necessary in order to explain the actively accreting close pairs, 
while the coexistence of close disk-less pairs with other accret- 
ing components in some systems argue in favor of no replenish- 
ment or a disrupted envelope. As for what is found for binaries, 
the disks seems to evolve in isolation at large separations and 
be affected by tidal truncation at small separations. 

Appendix A: Comments on individual objects 

RW Aur C : The suspected companion of RW Aur B (0712 
separation, K-band flux ratio of 0.024), detected in NIR by 
Ghez et al. ( 1993 1, was not detected. Our 5-sigma detection 
limit in the PSF-subtracted frames at the given position are 
0.03, 0.02 and 0.12 in Bry, H 2 and [Fell], respectively. This 
means that the companion would have been marginally de- 
tected in our data. As pointed out by Ghez et al. {1997b 1 and 
White & Ghez (2001 ) who did not detect the close pair in their 
HST optical images either, several reasons can explain our non- 
detection. First, orbital motion may have decreased the separa- 
tion below the resolution limit of NACO. The second expla- 
nation involves a drop of the flux of the companion since its 
detection in 1990 below our detection limit. Alternatively, the 



companion detected by Ghez et al. may in fact not be stellar 
but rather be the emission arising from shocked gas driven by 
a jet. RW Aur A is known to harbor an optical jet (Dougados 
et al. 2000 1 but no Herbig-Haro object is known to originate 
from RW Aur B. Another possibility, favored by White & Ghez 
J200T1 . is that this pair is spurious since the detection was ob- 
tained close to the detection limit and at relatively low SNR 
(Ghez et al. 1993 ). The reason of this false detection may how- 
ever be attributed to the presence of faint extended emission. 
In fact, the residuals of the PSF-subtraction clearly show some 
extended nebulosity around component B which could hardly 
be attributed to the difference of the AO-correction between the 
components of this l'/45 pair (Fig. IA.lt . 

Sz 19 B : The faint companion of Sz 19 A (K-band flux ra- 
tio of 0.023, Ghez et al. U997bl Chelli et al. [l995T l is not de- 
tected at the expected separation and PA (4.9±072, 202+3°). 
Our 5-sigma detection limits at that position are 0.05, 0.04 and 
0.04 in Bry, H2 and [Fell], respectively. Instead, we detect a 3- 
sigma emission at ~476 separation, PA~202°, which might be 
the companion. This is consistent with the values reported in 
RZ93 and the change of separation is within the uncertainties 
given by Ghez et al. dl997bi . 

HK Tau B : This companion, which was resolved as an edge-on 
disk (Stapelfeldt et al. 119981 Koresko[l998 ), is only marginally 
detected in our data at the 3-5 sigma level. Only the northwest 
nebula is detected at separations and PAs, as determined by 
the component centroids, 2'.'33, 171.6°, 2'.'29, 171.1° and 2731, 
171.1° in Bry, H 2 and [Fell], respectively. The position aver- 
aged over the three filters is 2731 ±0702 and PA=171.3±0.6°. 
The southeast nebula is a factor ~7 and ~8 fainter in integrated 
flux than its northwest counterpart in H and K, respectively 
(Koresko 1998 1, hence its non-detection. The position of the 
northwest nebula reported here is likely to be very similar to the 
position of the centroid of both nebulae, given their flux ratio. 
This position is slightly smaller than that reported in the visible 
HST/WFPC2 imaging of Stapelfeldt et al. (1T9981 . who mea- 
sured 2738, PA=172°. Since the flux ratios between nebulae in 
optical are similar to those in the near-infrared, the centroid in 
the optical would be coincident with that in the near-infrared. 
Leinert et al. ( 1993) also measured a somewhat larger separa- 
tion of 274±071 and PA=175±2°. In fact, we re-measured the 
position of the disk in the HST/WFPC2 images of Stapelfeldt et 
al. {1998 1 and find a better agreement with our derived NIR po- 
sitions with separations 2732±0701 and PA=171.7±0.2°. The 
derived uncertainties take into account the formal error be- 
tween the two images as well as plate scale and detector ori- 
entation uncertainties. We additionally confirmed our NACO 
separation and PA values with higher SNR NACO images in 
Ks which were available in the ESO archive and for which we 
derived a separation of 27318±07006 and PA=171.3±0.5°. All 
these values are consistent with an unchanged position of HK 
Tau B with respect to the primary. 

LkHa 263C: This edge-on disk (Chauvin et al. 120021 
Jayawardhana et al. 20021 was undetected for sensitivity rea- 
sons. Our 5-sigma detection limits at the expected position are 
0.06, 0.05 and 0.08 in Bry, H 2 and [Fell], respectively. These 
are about an order of magnitude larger than the actual flux ratio 



20 S . Correia et al. : A VLT/NACO Survey 




Fig.A.l. Result of the PSF-subtraction of RWAurB in Bry 
(top) and [Fell] (bottom). From left to right : primary and sec- 
ondary (1/3 power stretch display) and secondary minus pri- 
mary (linear stretch, upper cut is 10-sigma of the residual 
noise). FOV is l'.'4xl'.'4. 
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Fig.A.2. Astrometric measurements of SR24C. Offset posi- 
tions are shown with respect to SR 24 B. The expected position 
of a background star at the time of the NACO measurements 
(May 2004) is shown assuming it was located at the position 
measured by Simon et al. in August 1991. 

of the AC pair, determined by Jayawardhana et al. to be 0.007 
and 0.004 in H and K-band, respectively. 
Szl5 : Luhman (2004i, who recently produced a census of the 
Cha I region using also a comprehensive list of cross identifi- 
cations provided in Carpenter et al. ( 2002 ), concluded from the 
lack of Li absorption of Sz 15 (alias T19) and its low-extinction 
that it is actually a foreground star. Hartigan (1993) also noted 
the lack of Ha which precludes its identification as a YSO. 
PHa 30 AB : This system was indicated to be a background Be 
star based on an optical spectrum by Pettersson d!987i . 
SR 24N : By combining our new astrometric results with pre- 
vious measurements we can study the relative motion of this 
close pair. We include the 1991 measurements of Simon et 
al. ([1995 1 and measurements from 1999 unpublished archival 
HST/WFPC2 data in the FW606 and FW814 filters (HST pro- 
posal 7387, PI: K. Stapelfeldt). In addition, our measurements 
are confirmed by archival NACO data taken almost simultane- 
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Fig. A.3. Locus of the components of the LkHa 346 system in 
the (Teff, Mh) plane together with the PMS evolutionary tracks 
and isochrones of Baraffe et al. ( 1998 1. The components B and 
C are placed using their Mh and the assumption of coevality 
with either component A or D. 

ously to ours. The 1996 measurements of Costa et al. (|2000l 
are discarded from our analysis because of an obvious lack of 
accuracy. The locations of component C with respect to com- 
ponent B as a function of time are not consistent with a linear 
motion (Fig. lA.2| >. The NACO measurements unambiguously 
show that component C is in orbital motion around compo- 
nent B, rather than being an unrelated foreground/background 
object. Even if we don't consider the 1991 measurement, the 
resulting linear motion of velocity 19.0±6.5masyr~' at posi- 
tion angle 255.9°+13.1° is significantly different to the oppo- 
site of the proper motion of the system given by that of SR 24S 
which is 33.0+22.7 mas yr _1 at position angle 271.7°±39.3° 
(Ducourant et al. 2005 1. With the approach of the pericenter, 
it should be possible to derive a fairly accurate orbit in the 
next few years. Recently, McCabe et al. (2006 1 provide another 
measurement for an epoch intermediate between those of HST 
and NACO which fits quite nicely, although indicating that the 
marginally resolved HST measurement is of low quality. 
Sz 41 C : Known as CHX 20a in the optical identification of 
Einstein observations (Feigelson & Kriss 19891, it was later 
identified by Walter ( 1992 ) as a K0 III star, with radial velocity 
and Li 6707 A inconsistent with cloud membership. 
VW Cha : Like Brandeker et al. (200TJ. we did not see the 2'.'1 
separated candidate infrared companion found by Ghez et al. 
( 1997b I. Our 5-sigma flux ratio detection limits at that position 
are ~0.004 in all three narrow-band filters. 
LkHo' 346 : We give here details of the derivation of estimated 
spectral types and masses for the two new candidate compan- 
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ions (B and C) using their H-mag and the assumption of co- 
evality with the other components (A and D) of the system. 
For that purpose, we assume a distance of 160pc and assign 
A v =0.65 derived for LkHa 346 SE (Cohen & Kuhi [T979l here- 
after CK79) to components A, B, C and Ay=1.16 derived for 
LkHa 346 NW (CK79) to component D. The spectral type of 
component A is assumed to be M5, i.e. that of LkHa 346 SE 
(CK79), while the spectral type of component D is assumed 
to be K7, as the one of LkHa 346 NW (CK79). Converting to 
effective temperatures using the scale constructed by Sherry et 
al. J2004> gives T eff =3 125 K and 4060 K for components A and 
D, respectively. The masses and ages of these two components 
are derived by comparison with the Baraffe et al. ([1998 1 stel- 
lar evolutionary models in the Mh/Tej plane CFig. IA.3t . and 
are 0.28 + _° \ 6 2 M and < 1 Myr for component A, 0.8 M 
and 2 +j Myr for component D. Because the component A falls 
above the youngest isochrone, we constructed an extrapolated 
isochrone going through that locus from which the mass was 
derived. If we assume that the new candidate companions are 
coeval with component A, to which they are closest, this leads 
to a mass of 0.09 M and 0.09!°^ M for B and C, re- 
spectively. The derived range of T e ff corresponds to a spectral 
type of M6-7 for both B and C. If we assume coevality with 
component D rather than with component A we obtain masses 
and spectral types of ~0.2M Q and ~M4-6 for both compo- 
nents. 

Despite the large uncertainties involved in these mass and 
age estimates, a significant (> 1-sigma) non-coevality between 
component A and D is indicated. Alternatively, the relative lo- 
cus of component A with respect to component D could mean 
that the absolute magnitude is overestimated. This might be due 
in part to the H-band continuum excess which is not as negli- 
gable as the one in J-band (Greene & Meyer 1995). However, 
estimating the luminosity of the original two components of 
the binary (SE & NW, A and D component in our designation) 
via J magnitudes (or even visual magnitudes) and bolometric 
corrections leads to a similar age disparity. Using the luminosi- 
ties derived by CK79 increases the age of the SE component 
(~2Myr) but gives an unrealistically large age for the NW 
component of ~ 32 Myr. On the other hand, only the formal er- 
rors are included, and it might well be that model uncertainties 
could be largely responsible for the discrepancy. However, if 
one can consider the Baraffe et al. ( 1998 ) set of PMS tracks ac- 
curate enough, an overestimate of the luminosity of component 
A would have the consequence for the candidate companions B 
and C to be located in the substellar regime. Another potential 
source of error is the spectral type of component A. It might 
well be that the close proximity of components B and C have 
biased the spectral classification of CK79 to a later-type. This 
could then reconcile the age estimate of component A with that 
of component D. Accurate spectral types for all components of 
this system are required to further investigate this point. 
Hen 3-600 (TWA 3) : This binary exhibits a significant rela- 
tive motion of ~420mas at PA~143° since RZ93, which cor- 
responds to ~35.7masyr _1 . This relative motion is not con- 
sistent with the proper-motion measured for the system of 
137.8 mas yr 1 at PA=89.7° (Mamaiek l2005l and must there- 
fore originate from orbital motion. At the distance of the binary 



(50 pc, Tabled, the relative motion we measured corresponds 
to a linear relative motion of ~8.5 kms -1 , which is similar to 
a typical Keplerian velocity of 3.6kms _1 for such separation 
(~70 AU), considering a solar mass primary and a circular or- 
bit in the plane of the sky. In addition, including the 2000 mea- 
surement of Brandeker et al. ( 2003 B03 hereafter) shows that 
there is a significant change in the PA of the relative motion 
from 145° to 128°, from RZ93 to B03 and B03 to our mea- 
surement, respectively. We therefore conclude that the pair is 
physical and exhibits orbital motion. 

Appendix B: Edge-on disks in multiple systems 

Circumstellar disks in multiple T Tauri systems have only been 
imaged so far in the favorable detection case in which one disk 
is nearly seen edge-on. Three such systems exist: HKTauB, 
LkHa 262/263 C and HVTauC. The last two are in a system 
with more than 2 components (quadruple and triple, respec- 
tively). Disk inclinations are estimated to be ~85°, ~87° and 
-84°, while disk radius are -150 AU, -150 AU and -50AU, 
respectively. The following is an attempt to estimate lower and 
upper limits on the expected fraction of edge-on disks in sys- 
tems and compare the latter to the observed fraction in our sam- 
ple. 

If we consider a random distribution of disk inclinations 
then the probability for a disk to be at inclination ^ 85° is ~9%. 
Therefore, among all stars with a disk, 9% should have their 
disk oriented edge-on, if disks in multiple systems are sup- 
posed to be independently oriented. In the case in which disks 
of a pair are almost coplanar, this fraction should be increased 
to roughly 9 % of the pairs, neglecting systems with exactly 
coplanar disks whose detection in the edge-on orientation suf- 
fers from strong selection effects. Indeed, pairs with both disks 
edge-on would be difficult to detect as these sources would be 
faint, but if the disks are slightly misaligned (a few degrees 
from the edge-on orientation is sufficient for the star to become 
visible) these pairs would appear to contain a faint companion 
(like e.g. HKTauB) until high-angular resolution observations 
resolve its spatial extension. After removal of non-cloud mem- 
bers and probable chance projections our survey counts 42 dou- 
bles, 7 triples and 7 quadruples, i.e. a total of 133 components 
and 77 pairs. Two of the known edge-on disks are included in 
our sample (HKTauB and LkHa 262/263 C), which leads to 
a fraction of components and pairs of 2/133=1.5+1.1% and 
2/77=2.6+1.8%, respectively, i.e. more than 6- and 3-sigma 
lower than the above predicted fraction. 

However, the fraction of edge-on disks in our sample might 
be higher if we consider that : ( 1 ) they could not have been spa- 
tially resolved (this apply for the more distant objects and the 
smaller disk like the 50 AU radius HV Tau C), (2) the sensitiv- 
ity might not have been enough to detect the faintest one like 
LkHa 262/263 C which was not detected here (while HK Tau C 
has been, Sect.|Ajl, (3) the actual number of companions har- 
boring a disk in pairs is a fraction of the total number of pairs, 
since T Tauri pairs are known to have a significant fraction of 
mixed T Tauri types, i.e. disk-less companions (see Sect. l5.9> . 
On the other hand, although disks tend mostly to be coplanar 
with few exceptions as the latest polarimetric studies seems to 
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indicate (Monin et al. 2006a I, there are some indications for 
preferentially misaligned disks in high-order multiples (Monin 
et al. 2006b I. This would have the effect to reduce the predicted 
upper-limit on the fraction of edge-on disks derived above as- 
suming coplanar disks for all pairs, irrespective of their multi- 
ple system memberships ( the 9 % of the pairs). While the exis- 
tence of a deficit of edge-on disk in our sample is difficult to es- 
tablish firmly, the fact that 2/3 of the known edge-on disks are 
found in high-order multiple systems suggest that misaligned 
disks are more likely to occur in these systems compared to bi- 
naries because of internal dynamical evolution in unstable mul- 
tiple systems (Monin et al. 2006b). 

It is also worth to remark that an infrared companion (IRC, 
e.g. Koresko et al. l 19971 1 candidate, namely FVTauD, is in- 
cluded in our sample of high-order multiple systems. This is 
consistent with the frequency of ~ 10% with which this class of 
objects is thought to occur in the T Tauri multiple stars popula- 
tion. Together with V 773 Tau D, WL 20 S and T Tau Sa, there 
are an increasing number of IRCs 7 that are members of high- 
order multiple systems, which may possibly hint at a relation 
between the IRC phenomenon and the degree of multiplicity. 
This might favor, at least for some of them, the scenario pro- 
posed some years ago by Reipurth (2000 1. 
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